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ABSTRACT 

During the embedded phase of pre-main sequence stellar evolution, a disk forms from the dense envelope while an accretion-driven 
outflow carves out a cavity within the envelope. Highly excited (E' = 1000-3000 K) H2O emission in spatially unresolved Spitzer/IRS 
spectra of a low-mass Class object, NGC 1333 IRAS 4B, has previously been attributed to the envelope-disk accretion shock. 
However, the highly excited H2O emission could instead be produced in an outflow. As part of the survey of low-mass sources in 
the Water in Star Forming Regions with Herschel (WISH-LM) program, we used HerschelfPACS to obtain a far-IR spectrum and 
several Nyquist-sampled spectral images to determine the origin of excited H2O emission from NGC 1333 IRAS 4B. The spectrum 
has high signal-to-noise in a rich forest of H2O, CO, and OH lines, providing a near-complete census of far-lR molecular emission 
from a Class protostar The excitation diagrams for the three molecules all require fits with two excitation temperatures, indicating 
the presence of two physical components. The highly excited component of H2O emission is characterized by subthermal excitation 
of ~ 1500 K gas with a density of ~ 3 x 10* cm"', conditions that also reproduce the mid-IR H2O emission detected by Spitzer. On 
the other hand, a high density, low temperature gas can reproduce the H2O spectrum observed by Spitzer but underpredicts the H2O 
lines seen by Herschel. Nyquist-sampled spectral maps of several lines show two spatial components of H2O emission, one centered 
at ~ 5" (1200 AU) south of the central source at the position of the blueshifted outflow lobe and a second centered on-source. The 
redshifted outflow lobe is likely completely obscured, even in the far-IR, by the optically thick envelope. Both spatial components of 
the far-IR H2O emission are consistent with emission from the outflow. In the blueshifted outflow lobe over 90% of the gas-phase O is 
molecular, with H2O twice as abundant than CO and 10 times more abundant than OH. The gas cooling from the IRAS 4B envelope 
cavity walls is dominated by far-IR H2O emission, in contrast to stronger [O I] and CO cooling from more evolved protostars. The 
high H2O luminosity may indicate that the shock-heated outflow is shielded from UV radiation produced by the star and at the bow 
shock. 

Key words. Stars: protostars; ISM: jets and outflows; Protoplanetary disks; Infrared: stars; Techniques: spectroscopic 



^ . 1. INTRODUCTION heated by shocks and energetic radia tion from the central 

rrj star ( e.g.lSneU et all [19 80: Scaans et aTlll995HArce & SargentL 

^ ^ During the embedded phase of pre-mam sequence stellar evo- |^|van Kemp^ITetTl., 2009,; JoMnetaLillOH*. On theothTr 



H , lution, the protoplanetary disk forrmout of a dense molec- j^and, observational evidence for the disk-envelope interactions 

. ^ . ular envelope (e.g. [ Terebeyet a^ [19841 [ Adams et a^ [1987[). ^as been SDai-se. lVelusamv eTdl ([S) detected methanol emis- 

Meanwhde, as the protostar builds up most of its mass, it ^^^^ ^^.^^ L1157 on scales of ~ 1000 AU, spatially-extended 

drives powerful, coUimated outflows into the dense envelope beyond the point-hke continuum emission, and argued that 

(e.g. [ BontempsetaH \1996). These processes together eventu- j^e kinematics suggest that the emission is produced at the 

ally cause the envelope to dissipate and set the initial conditions disk/envelope interface. Watson et al. ( 2007) detected emission 

for disk evolution and planet formation. j„ highly excited (E' = 1000 - 3000 K) H2O lines from the NGC 

At the interfaces between the outflow and envelope and be- 1333 iraS 4B system and attributed the heating to an envelope- 

tween envelope and disk, shocks can heat the gas and poten- disk accretion shock within 100 AU of the star. These obser- 

tially produce detectable emission. The well-studied outflow- yations off^er two dilferent interpretations for disk-envelope in- 

envelope interactions produce an outflow cavity with walls teractions, w ith material either entering the disk on large scales 

(IVisser et all 2009; Vorobvov, 2010) or raining onto the disk 

* Herschel is an ES A space observatory with science instruments pro- at small radii (Whitney & Hartmann. ,1993 ). However, a persis- 

vided by European-led Principal Investigator consortia and with impor- jgnt complication in interpreting spatially and spectrally unre- 
tant participation from NASA. 
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Fig. 1. The location of the 5x5 spaxel array (purple spectral map) 
again st an 4.5 jum image obta ined with SpitzerflRAC (grayscale) 
from 'j0rgens en et"an (l200d) . The spectral map shows that the 
continuum-subtracted emission in the 0-H2O 616 - 5o5 82.03 fim 
(E' - 643 K) line is produced mostly in the blueshifted out- 
flow lobe. The sub-mm positions of IRAS 4B and IRAS 4A are 
marked with blue asterisks. 



solved emission as coming from a compact disk-like structure 
is that outflows can also produce bright emission in highly ex- 
cited lines. Molecular emission is a dominant coolant of outflow- 
envelope interactions, and H2O emission is particularly sensitive 
to shocks in outflows (e.g. [Nisini et al., 2002; van Kempen et al., 

IMM lEiMetll lloio). 

The NGC 1333 IRAS 4B system (hereafte r IRAS 4B) is a 
Class YSO (d = 235 pc. lHirota eta"ni2008h with a 0.24 Mg 
disk that is deeply embedded (Ay ~ 1000 mag) within a 2.9 
M0 envelope (J0rgensen et al. 2002, 2009). Compact outflow 
emission is detected in many sub-mil l imeter (sub-mm) molec- 
ular Unes dPi Francesco et all I2OOII: lJ0rgensen et al.L l2007h . 
Near-IR emission in all four Spitzer/IRAC bands (3.5, 4.5, 
5.8, and 8.0 /jm) is located in the blueshifted outflow lobe 
(J0rgensen et al. 2006, see also Choi et al. 2011), offset by 
~ 6" south of the peak of interferometric sub-mm continuum 
emis sion. The redshif ted outflow is seen in sub-mm line emis- 
sion (|j0rgensen et al.L [2007) but is not detected in the near- or 
mid-IR because the outflow is located behind IRAS 4B and hid- 
den by the high extinction of the envelope. In the near-IR, the 
Spitzer/IRAC photometry is dominated by H2 emission (Arnold 
et al. 201 1; see also Neufeld et al. 2008), with some contribution 
of CO fundamental emission to the 4.5 yum bandpass (Tappe et 
al. 2011; see also Herczeg et al. 2011). Excited water emission 
from the NGC 1333 IRAS 4 system, including both IRAS 4A 
and 4B, was detected with ISOfLWS, but with too low spatial 
resolution to attribute the emission to any component in the sys- 
tem (Giannini et al., 2001). H2O maser emission has also been 
seen from dense gas associated with the IRAS 4B outflow, al- 
though with a different position angle than the molecular outflow 
(Rodriguez et al., 20021: iFuruva et al.L 120031; [Marvel et aUl2008l; 
Ibesmurs et al.1. 1200 9). 

In a sampl e of low-reso l ution Spitzer-flRS spectra of 30 
Class objects. IWatson et"an ( l2007h found water emission from 
highly-excited levels in only IRAS 4B. The lines have upper 



levels with high excitations (1000 - 3000 K) and high critical 
densities (~ 10^' cm""^). Watson et al. inferred that the emitting 
gas has a high density and argued that the high density indi- 
cates that the emission is produced in a ~ 2 km s"' accretion 
shock at the envelope-disk interface. They also argued that the 
emission was detected only from IRAS 4B because the view- 
ing angle may be well-aligned with the outflow, allowing a clear 
view of the embedded disk, and because the timescale for such 
high envelope-disk accretion rates may be short. H2O emission 
has since been detected in at least two other components of the 
IRAS 4B system: (1) narrow (~ 1 km s"', with rotation signa- 
tures) p-Hj^^O 3 13 - 220 (E' - 204 K) emission in a spatially 
compact region, likely a (pseudo )-disk of ~ 25 AU in radius 
( lJ0rgensen & van DishoecM.l20I0h ; and (2) broad (FWHM- 24 
km s~'), spatially unresolved emission from low-excitation H2O 
(£■' - 50 - 250 K) lines, wh ich are consistent with an outflow 
origin ( Kristensen et al.L 1201 0) but too broad for the ~ 2 km s ' 
velocity expecte d of envelope gas in free-fall striking the disk at 
25-100 AU (e.g. lShuetal.lll977h . 



The highly-excited mid-IR emission and the broad line pro- 
files of lower-excitation lines could be reconciled if either (a) 
the high- and low-excitation H2O line emission originates in dif- 
ferent locations, if (b) models for the envelope-disk accretion 
shock underpredict line widths, or if (c) both the high- and low- 
excitation H2O lines are produced in the outflow. In this paper, 
we analyze a HerschelfPACS far-IR spectral survey and spectral 
imaging of IRAS 4B to resolve the discrepancy in the differ- 
ent possible origins of H2O emission from IRAS 4B. The far-IR 
spectrum of IRAS 4B is as rich in lines as its mid-IR Spitzer/IRS 
spectrum. In Nyquist-sampled maps, the H2O emission is spa- 
tially offset from the continuum emission in the direction of 
the blueshifted outflow. The excitation of warm H2O lines de- 
tected with Herschel/PACS and with SpitzerflRS together can 
be explained by emission from an isothermal slab. An additional 
physical component(s) is located on-source and likely traces ma- 
terial closer to the base of the outflow. No evidence is seen for 
an envelope-dis k accretion s hock. This result is consistent with 
a conclusion bv lTappe et al.l l2012l) . obtained contemporaneous 
to the results in this paper, that the H2O emission seen in the 
Spitzer/IRS spectra also coincide with the southern outflow po- 
sition. We discuss the implications of these results for outflows 
and for the prospects of observationally studying disk formation 
with H2O lines. 



2. OBSERVATIONS AND DATA REDUCTION 

We obtained far -IR spectra of NGC 133 3 IRAS 4B (03''29"'12 '.0 
-h31°13'08'.'1; lJ0rgensen et all IIoOTI) on 15 -16 Marc h 2011 
with the R^CS instru ment on Herschel (Pilb ratt et all I2OIOI: 
Poglitsch et a ll 1201 ()h as part of the WISH Key Program 
(van Dishoeck et al. □20111) . The observations presented here 
consist of a complete scan of the 52-208 /im spectral range 
and deep, Nyquist-sampled spectral maps in four narrow (~ 
0.5 - 1 jjm) wavelength regions. Each PACS spectrum in- 
cludes observations of two different nod positions located at 
+3' from the science observation to subtract the instrumental 
background. All PACS data were reduced with HIPEv6.1 Jotj, 
. We supplemented the PACS observations with re-reduced 



'2010) 

archival Spitzer/IRS spectra that were previously analyzed by 
IWatson et al] ( l2007h . Details of the observations and reduction 
are described in the following subsections. 
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Fig. 2. The FWHM from initial fits to strong lines in our full scan 
of the 52-208 jum wavelength range. The lines are spectrally un- 
resolved and provide a measure of the spectral resolution. The 
line widths for our final fits were set by linear fits to the FWHM 
at>100 and<100jum. 



2.1. Complete Herschel/PACS far-IR spectrum 

The 52-208 fim spectrum of IRAS 4B was obtained in 2.7 hrs 
of integration with PACS. PACS observed IRAS 4B simultane- 
ously in the first order > 100 /im and in the second order at < 100 
yum. The grating resolution varies between R = 1000 - 2000 at 
> 100 jum and R = 3000 - 4000 at < 100 fim. Large grating 
steps are used, so the spectrum is binned to ~ 2 pixels per reso- 
lution element. A spectral flatfielding was applied to the data to 
improve the S/N. The observed flux was normalized to the tele- 
scopic background and subsequently calibrated from observa- 
tions of Neptune, which is used as a spectral standard. Appendix 
A describes a new method to calibrate the flux in PACS spectra 
between 97-103 fim and above > 190 fim. The relative flux cal- 
ibration is accurate to ~ 20% across most of the spectrum. 

The spectral scan produced a single 5x5 spectral map over 
a 47" X 47" field-of-view. The central spaxel is centered at the 
location of the sub-mm continuum peak. An adjacent spaxel lo- 
cated 9'.'4 to the SE (PA=249°) is centered on the blueshifted 
outflow position (Fig.[TJ. Most of the line and continuum flux is 
located in these two spaxels. The southern outflow of NGC 1333 
IRAS 4 A is located in the NW corner of the array. 

Extracting line fluxes requires an assessment of the distri- 
bution of flux on the detector caused by both the point-spread 
function of Herschel and the spatial extent of the emission. For 
a well-centered point source, the encircled energy of in a single 
spaxel is ~70% at < 100 fim wavelengths and declines to ~ 40% 
near 200 yum. However, the signal-to-noise decreases if the spec- 
trum is extracted from many spaxels. The emission line spec- 
trum is obtained by adding the flux from the central spaxel and 
the outflow spaxel. The line fluxes are subequently corrected for 
light leakage and the spatial extent in the emission by compar- 
ing line fluxes from this 2-spaxel extraction with the line fluxes 
extracted from a 3 X 3 spaxel area centered on the central spaxel. 
The line flux ratio between the 2-spaxel and 3x3 spaxel extrac- 
tion was calculated for strong lines of all detected molecules. 
The fluxes from the 2-spaxel extraction are then divided by a 
wavelength-dependent correction that is 0.76 at < 100 fim and 




20 10 -10 -20 

ARA (arcsec) 



Fig. 3. The 63 fim continuum emission (colors) and the SCUBA 
450 /im emission map (contours), with the location of IRAS 4A 
and IRAS 4B marked as white asterisks. The PACS maps are 
shifted so that the centroid of the 190 fim continuum emission is 
located at the position of IRAS 4B obtained from sub-mm con- 
tinuum interferometry. The SCUBA map is shifted to the same 
position. 

then decreases linearly to 0.6 at 180 jim. The wavelength de- 
pendence of this correction includes both the point-spread func- 
tion of Herschel and the wavelength dependence in the spatial 
distribution of the detected emission. Finally, hi gh signal-to- 
noise PACS spectra of the point source HD 100546 (IS turm et all 
2010) were then used to correct for emission leaked beyond the 
3x3 spaxel area. This approach assumes a similar spatial dis- 
tribution for all molecules and for lines at nearby wavelengths, 
and introduces a ~ 10% uncertainty in relative fluxes. The over- 
all uncertainty in flux calibration is ~ 30%. The typical noise 
level in the two-spaxel extraction (~ 9.4 x 18.8") is ~ 0.4 Jy per 
resolution element in the extracted continuum spectrum. 

All lines in the complete spectral scan are spectrally un- 
resolved. The lines were initially fit with Gaussian profiles, 
with the central wavelength, line width, and amplitude and 
a first-order continuum as free parameters. The instrumental 
line widths were then calculated from first-order fits to the 
wavelength-dependent spectral resolution, obtained from strong 
lines (Fig.|2]i. Our final fluxes are obtained from fitting Gaussian 
profiles to each line, with the line width set by the calculated in- 
strumental resolution at the given wavelength. For strong lines, 
the median centroid velocity is 30 km s ' with a standard devi- 
ation of 24 km s"' . The absolute wavelength calibration is accu- 
rate to ~ 50 km s"' and is limited by the spatial distribution of 
emission within each spaxel. 

2.2. Herschel/PACS Nyquist-sampled spectral Imaging 

Nyquist-sampled spectral maps of narrow spectral regions at 
54.5, 63.3, 108.5, and 190 //m were obtained in a total of 2.2 
hr of integration time. These maps were obtained from a 3 x 3 
raster scan with 3" steps, yielding spatial resolutions of ~ 5", 5", 
8", and 10", respectively. Small grating steps were used to fully 
sample the spectral resolution. The final spectrum is rebinned 
onto a wavelength grid with ~ 4 pixels per resolution element. 
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Fig. 4. The combined HerschelfPACS and Spitzer/IRS continuum-subtracted spectra of IRAS 4B, with bright emission in many 
H2O (blue marks), CO (red marks), OH (green marks), and atomic or ionized lines (purple). The Spitzer spectrum is multiplied 
by a factor of 10 so that the lines are strong enough to be seen on the plot. The inset shows the combined spectrum including the 
continuum. 



The spectral maps include CO, H2O, and [O I] lines listed 
in Table [T] In a 3 x 3" area, the typical rms is ~ 0.02 Jy per 
resolution element at 63 fxm and ~ 0.01 Jy per resolution ele- 
ment at 108 yum. This sensitivity level is better than that from the 
full spectral scan because of different integration areas and much 
longer integration times in each resolution element. 

The data cubes from the Nyquist-sampled maps were repro- 
jected onto a normal grid of right ascension and declination. In 
the automated calibration the 190 yum continuum emission is off- 
set from the o bject position', as measured from sub-mm interfer- 
ometry ( J0rg ensen et all l2007h . Each map is shifted in position 
by 3'.'3 E and 1'.'5 N so that the location of the 190 jum con- 
tinuum from our observations matches the peak location of the 
sub-mm emission. After the shift, the 63 and 108 /vm contin- 
uum emission from IRAS 4B and IRAS 4A (located in the NW 
edge of the map) are well aligned with the Sp/fzer-MIPS 70 yum 
emission. No significant offset was measured in the staring ob- 
servation of the full PACS SED, which was obtained with a new 
pointing one day after the Nyquist-sampled map. Because the 
same spatial shift is applied to both the line and far-IR contin- 
uum emission maps and because the shift moved both the line 
and continuum closer to the sub-mm continuum peak, the result 
that the line emission is spatially offset from the sub-mm contin- 
uum is robust to the pointing uncertainty. 

The southern nod position in the map has spatially-extended 
[O I] emission in the southwest portion of the map. The [O I] 



' The optical guider observations onboard Herschel are typically ac- 
curate to ~ 1". However, observations of nearby high-extinction re- 
gions, including NGC 1333 IRAS 4B, have few optical guide stars, 
which may introduce pointing offsets when the few guide stars are dis- 
tributed asymmetrically in the field-of-view. 



emission is therefore measured from only the northern nod po- 
sition. Inspection of 5/7/fzer/MIPS 70 //m maps at the two nod 
positions does not indicate the presence of a strong, point-source 
continuum emission that could otherwise corrupt the continuum 
map for IRAS 4B. The 54 /im continuum map has low S/N and 
is not used. 

The two H2O and [O I] lines at 63 yum have spectral widths 
of 1 10 km s which places an upper limit of 65 km s ' on the 
intrinsic line width. That the line widths are broader than the in- 
strumental resolution of ~ 3300 is not significant because emis- 
sion that is spatially extended in the cross-dispersion direction 
can broaden the spectral line profile, as with any other slit spec- 
trograph. 

2.3. Spitzer/IRS Spectrum 

The Spitzer/IRS spectra of IRAS 4B were originally presented 
by Watson et al. (2007). We re-reduced the spectrum following 
the procedure described by Pontoppidan et al.l (l2OI0l) . The LH 
slit width is 1 1". The flux extraction region of 5-10" across the 
wavelength region was selected to o ptimize the fi nal signal-to- 
noise in the limit of a point source ( Hornel 1198 6*). The Spitzer 
and Herschel observations cover similar regions on the sky. The 
relative flux calibration between Spitzer and Herschel spectra is 
likely uncertain by ~ 30% flux. 

3. RESULTS 

3.1. Far-IR Spectrum of IRAS 4B 

The far-IR PACS spectrum of IRAS 4B is the richest far-IR spec- 
trum of a YSO to date. A forest of high signal-to-noise CO, 
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Fig. 6. Left and Middle: Contour maps of H2O (yellow), [O I] (red, left), and CO (red, middle) emission compared with continuum 
emission (color contours) at 63.3 fim (left) and 108.5 fim (middle). The insets show the spatial extent of continuum and line emission 
in the different components along the N-S outflow axis. Right: a comparison of the location of emission in H?0 63.32 um (b lue), 
108.5 fim continuum (red), Spitzer-IRAC 4.5 jum photometry (color contours: 'J 0rgensen et al.L I2OO6I : iGutermuth et al.L 12008 ). All 
contours have levels of 0.1,0.2,0.4,0.8, and 1.6 times the peak flux near IRAS 4B. The asterisks show the sub-mm positions of IRAS 
4B near the center of the field and IRAS 4A in the NW corner Some noise in these maps is suppressed at empty locations far from 
IRAS 4A and IRAS 4B. 



: Continuum 

:/ 


spa 
Ou 


xel 1 jt_Lujj_ 

$ 

tflow spaxel 














il 
k 






L 




t** - 
i>f\ - 



50 100 200 

Wavelength (^m) 



Fig. 5. Top: Spectra extracted separately from a spaxel centered 
on the sub-mm continuum (red) and a spaxel offset by 9'.'4 to the 
S and centered on the outflow position (blue). The outflow po- 
sition is dominated by line emission while the central object is 
dominated by continuum emission, indicating that the lines and 
continuum are spatially offset. Bottom: The same two spectra 
after continuum subtraction show that the on-source line emis- 
sion gets much weaker to short wavelengths. 



H2O, and OH lines that provide a full census of far-IR molec- 
ular emission that can be detected from low-mass YSOs (Fig. 4; 
see also Figs. lD.l]|D.2l in the Appendix). Lines were discovered 



in the spectrum following a biased search for emission at the 
wavelengths for transitions of common species and an unbiased 
search for narrow features that peak above the noise level. 

A total of 115 distinct emission lines are detected and iden- 
tified from IRAS 4B (Table 3). All strong lines are identified. 
Several tentative detections of weak lines are unidentified and 
discussed in Appendix B. No Hj^O or '^CO emission is detected 
in the PACS observations, with typical flux limits in the strongest 
expected lines of ~ 0.03 times the observed flux of the main iso- 
topologue. The [O I] 145.5 jum line is not detected, with a 2cr 
flux Hmit 1.2 X 10-21 W cm ^. Lines of OH+ and CH+, HD 56 
and 1 1 2 //m, [N II] 1 2 1 . 8 and 205 .2 //m, [C II] 1 57 .7 fim, and [O 
III] 88.7 fj.m are also not detected. 

Figure[T]shows a spectral map of the 0-H2O 616 - 5o5 82.03 
fim (E' = 643 K) line emission overplotted on a SpitzerflRAC 
4.5 jum image of IRAS 4B (and IRAS 4A). The line emission is 
located at the position of the near-IR emission in the blueshifted 
outflow lobe, south of the central source of IRAS 4B. On the 
other hand, most of the continuum emission is located in the cen- 
tral spaxel, consistent with the location of the sub-mm contin- 
uum emission. Figure |5] demonstrates that the equivalent width 
of far-IR lines is much larger in the outflow spaxel than in the 
continuum spaxel, indicating a spatial offset between the line 
and continuum emission. The line emission in the central spaxel 
mostly disappears at < 70 /um. All lines in the R4CS spectrum 
are spatially offset from the continuum emission. The line fluxes 
are measured based on the summation of these two bright spax- 
els and correction for spatial extent (see §2.1). 

3.2. PACS Mapping of H2O Emission from IRAS 4B 

Figure |6] and Table [1] compare the spatial distribution of the 
continuum flux with H2O, CO, and [O I] line fluxes obtained 
from the Nyquist-sampled spectral maps. The far-IR continuum 
emission is centered on-source, at the location of the sub-mm 
continuum, while line emission is centered to the south in the 
blueshifted outflow. Figure |7] shows a cartoon version of the ap- 
proximate location of line and continuum emission and the mor- 
phology of IRAS 4B. In the following analysis, we simplify the 
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Fig. 7. A cartoon showing the location of different emission 
components from IRAS 4B, based in part on Fig.|6l 



analysis by assuming that emission consists of two unresolved 
point sources, one at the outflow position and one at the sub- 
mm continuum peak, and are subsequently fit with 2D Gaussian 
profiles. More complicated spatial distributions would be unre- 
solved in our maps. 

The right panel of Figure ^demonstrates that the location of 
the warm H2O coincides with the Spitzer/IRAC 4.5 /vm imag- 
ing. The two H2O lines near 63.4 fim, 0-H2O 818 - Iqj and p- 
H2O 808 - 7,7, (£" = 1070 K) are centered at 5.2 + 0.2" from 
the 63 jum continuum and are spatially extended relative to the 
continuum emission (assumed to be unresolved for simplicity) 
by FWHM= 6.1 + 1.0". About 70% of the emission is produced 
at the southern outflow position (Fig. [8]). The flux ratios for the 
two H2O 63.4 jum lines are similar at both the on-source and 
off-source positions (Fig.|9l). In the 108.5 yum spectral map, both 
the 0-H2O 221 - lio (£' = 194 K) and CO 24-23 (£' = 1524 
K) emission are centered 2.5 ± 0.4" south of the 108.5 jum con- 
tinuum emission and are spatially-extended in the outflow di- 
rection by 13.6 + 0.7", relative to the extent of the continuum 
emission. The larger spatial extent and smaller offset in the 108 
fj.m lines both indicate that the outflow component contributes 
~ 40% of the measured line flux. The spatial differences may be 
interpreted as differential extinction across the emission region, 
discussed in the next subsection. 

The 54 jum maps are noisy because PACS has poor sensitivity 
at < 60 fim. The 0-H2O 532 - 5o5 54.507 fim {E' = 732 K) 
emission is offset by 5.9 + 0.4" south from the the peak of the 
sub-mm continuum emission. The CO 49-48 53.9 fim emission 
(£" = 6457 K) is offset 2.9 ± 0.4" south, between the peak of the 
sub-mm emission and the bright outflow location. The highly- 
excited CO emission is produced in a different location than the 
highly-excited H2O emission. 

The [O I] emission is offset by 3.7 + 0.3" at PA=168°, just 
west of the outflow, and is spatially extended by ~ 7'.' 1 ± 1 .0. 

3.3. Extinction Estimates to tine Centrai Source and Outfiow 
Position 

The extinction to different physical structures within the IRAS 
4B system depends on how much envelope material is present in 
our line of sight. In this subsection, we discuss how these differ- 
ent extinctions affects the emissi on that is seen. The ext i nction 
law used here is obtained from iWeingartner & Draind (1200 lb 
with a total-to-selective extinction parameter Rv-5.5, typical of 



dense regions in mole cular clouds ("e.g. llndebetouw et all 120051 : 
IChap man et ali l2b09l) . Appendix D includes a discussion of how 
extinctions may affect the molecular excitation diagrams. 

The strength of the near-IR emission in the outflow 
(|j0rgensen et all l2006l) indicates that extinction must be low to 
at least some of the outflow position. Spherical models of the 
dust continuum indicate that the central protostar is surrounded 
by Ay - 1000 mag. (J0rgensen et al., 2002), so any emission 
from the redshifted outflow lobe may suffer from as much as 
Ay ~ 2000 mag. of extinction. Any additional extinction would 
have likely introdu ced asymmetries i n the H 2O line profiles that 
were presented in Kristens en et alj (|2010|) . Depending on the 
wavelength and spatial location of the emission, the far-IR emis- 
sion line fluxes may be severely affected by extinction. 

The H2O 54.5 and 63.4 jjm lines have a different spatial dis- 
tribution than the H2O 108.1 //m line. If we assume that the on- 
source and off-source emission both have similar physical con- 
ditions, then the ratio of the H2O 63.4 to 108.1 /urn line lumi- 
nosities should not change with position. In this scenario, the 
different locations for the detected flux is caused by differential 
extinction across the emission area. An average extinction to the 
on-source H2O component of Ay ~ 700 mag. would reduce the 
fractional contributions from the on-source and off-source loca- 
tions observed values. This extinction may be the combination 
of a lightly-extincted region on the front side of the protostar 
and a heavily-extincted region on the back side of the protostar. 

An independent estimate of the extinction can also be made 
from the flux ratio of [O I] 63.18 to 1 45.5 nm lines, which 
is typically observ ed to be about 10 dGiannini et al 
iLiseauet all 12006*). The undetected [O I] 145.5 fim line flux is 
less flian 10% of the 1.8 x 10"^° W cm ^ flux in the [O I] 63.18 
fim line. If we conservatively assume that the true ratio is 30, 
then we estimate Ay < 200 mag. to the [O I] emission region. 
Some additional [O I] emission could only be hidden behind a 
high enough extinction (Ay ~ 4000 mag.) to attenuate emission 
in both the 63.18 and 145.5 jim lines. Therefore, the effect of 
extinction on the [O I] luminosity is likely not too significant for 
the southern outflow lobe, where [O I] emission is seen. 

3.4. Spitzer/IRS Spectrum and Broadband Images of IRAS 
4B 

The Spitzer/IRS spectrum of IRAS 4B includes emission in lines 
of highly-excited H2O and OH, plus [ S I] and [Si II ] Our Une 
identification mostly agrees with that of lWatsonetaP (l2007h for 
H2O lines, with some modifications to account for the identifi- 
cations of OH lines (Table |2] and Figure [10] see also Tappe et 
al. 2008). The H2O line identification was informed from line 
intensities predicted by RADEX modelling (see §4.1) of both 
th e low-density case d iscussed here and the high density-case 
of IWatson et"al] (l2007h . All lines with significant detections are 
identified. 

An analysis of the SpitzerfMIPS 24 jjm image, with sensitiv- 
ity from 20-3 1 fim, helps us place limits on the amount of H2O 
emission that could arise on-source. Convolving the filter trans- 
mission curve with the Spitzer/IRS spectrum from Watso n et al.l 
(2007) indicates that ~ 24% of the Ught in the MIPS 24 //m 
bandpass is in molecular emission (mostly H2O), 17% in the [S 
I] 2 5.24 yum line, and 59% in th e continuum. 

|j0rgensen & van Dishoeckl (|2010) found that the emission 
from IRAS 4B in the SpitzerfMIPS 24 fim images is offset 
from the pe ak emission of the sub-mm continuum (see also 
IChoi&Led bOll). We measure that the emission is centered 
at 5'.'2 S and 0'.'9 E from the central source, consistent with the 
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Fig. 9. The H2O 63.4 fim spectral region extracted from the on- 
source position (red) and from the blueshifted outflow position 
(blue), and over the entire spectral map (black). No significant 
differences are detected in the ratio of the two lines, which indi- 
cate that the two lines are optically thin at both locations. 



Fig. 8. Top: A cross cut of the flux at several wavelengths along 
the N-S outflow axis. The 63 //m continuum is centered at the 
peak of the sub-mm continuum while the 4.5 fim and 24 fim 
photometry are centered in the outflow. The H2O 63.3 fj.m lines 
are produced primarily at the outflow location while the H2O 
108.1 fim line has similar on- and off-source contributions. The 
(0,0) position is de fined here by the peak o f the sub-mm contin- 
uum emission from |j0rgensen et aTl (l2007h . Bottom: The spatial 
cross cut of H2O 63.3 jum and MIPS 24 /im emission, shown 
as the combination of two unresolved Gaussian profiles located 
at the on-source position (red dashed line) and the blueshifted 
outflow position (blue dashed line). 



location of the outflow emission, and is spatially extended by 
~ 5'.'5 in the north-south direction, along the outflow axis. An 
additional component is present at the location of the sub-mm 
continuum peak. The Spitzer/MlPS emission is assumed here to 
be a combination of emission from two unresolved sources, one 
at the sub-mm continuum peak and one at the position of the 
Spitzer IRAC 4.5 yum emission located 6'.'2 S and 0'.'4 E. From 
fitting two dimensional Gaussian profiles to the image, the com- 
ponent at the blueshifted outflow lobe accounts for 78% of the 
Spitzer/MlPS 24 fim emission and the sub-mm point source ac- 
counts for the remaining 22% of the emission (see Fig.[8]for the 
fit to the image collapsed onto the outflow direction). Some addi- 
tional Spitzer/MlPS emission is located at 20" S of the sub-mm 
continuum peak and is ignored here. 

In the Spitzer/IRAC images of emission between 3.8-8 fim, 
the emission is located entirely at the outflow position. In con- 
trast, the 63 pm continuum emission is located mostly on the 
central source at the sub-mm continuum position. Much of the 
20-3 1 yum continuum emission must be located at the outflow po- 
sition, but some continuum emission could also be located on the 
central source. Given the fraction of emission located on-source 
(22%) and the relative contributions of molecular lines (24%) 
and continuum to the SpitzerfMlPS photometry, the MIPS map 
could be consistent with an on-source location of H2O emission 
only if the continuum emission is located entirely at the outflow 
position. 



4. EXCITATION OF MOLECULAR EMISSION FROM 
IRAS 4B 

The spatial distribution of the highly excited H2O emission in 
the PACS observations places the bulk of the highly excited far- 
IR emission at the outflow position. In this sections, we ana- 
lyze the excitation of the H2O lines in detail to demonstrate that 
the highly excited H2O emission in both the Herschel/PACS and 
Spitzer/IRS spectra can be explained with emission from a sin- 
gle isothermal, plane-parallel slab. We subsequently analyze CO, 
OH, and [O I] emission from IRAS 4B. Although the H2O emis- 
sion region is likely complicated and includes multiple spatial 
and excitation temperature components, our simplified approach 
is able to reproduce the highly excited H2O lines. The properties 
of this slab are the combination of the spatially offset outflow 
component and the on-source component. We lack sufficient spa- 
tial resolution throughout most of the spectrum to analyze the 
excitation of the two components separately. 

Figure [TT] shows excitation diagrams for H2O, CO, and 
OH emission^. Without considering sub-thermal excitation, each 
molecule requires two temperature components to reproduce the 
measured fluxes. For convenience, the two components for each 
molecule are called "warm" and "cool", however this terminol- 
ogy applies separately to each molecule^. The warmer compo- 
nent of CO may not be related to the warmer component of H2O 
or OH. For the temperature and density derived below, the two 
apparent excitation temperatures for OH and H2O could even 
be produced by a single component, with the warm and cool 
regimes resulting from subthermal excitation. Table [3]describes 
the excitation temperatures, molecular column density, and lu- 
minosity for fits to these diagrams. In the following subsections 
we discuss the excitation of the H2O in detail, and briefly de- 
scribe the excitation of CO, OH, and O. RADEX models of H2O 
are used to fit only the higher excitation H2O lines because the 
lower excitation lines are optically thick and difficult to use to 
infer physical conditions of the emitting gas. The emitting area, 
temperature, and density derived from the fits to the observed 
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Fig. 10. OH emission lines (vertical dashed lines) in the 
Spitzer/IRS spectrum of IRAS 4B. The lines expected to be 
strongest are either detected (Table |2]l or are blended with lines 
of other species. 



H2O lines are assumed to also apply to CO, OH, and O for sim- 
plicity. 

4.1. RADEX Models of H2O Emission 

The H2O line emission extends to high energy levels, with level 
populations indicating an excitation temperature of 220 K but 
with significant scatter Among the highly excited lines, the 
H2O excitation diagram does not show any break at high en- 
ergies that would indicate multiple excitation components. The 
highly-excited H2O levels have high critical densities (~ 10" 
cm"^). The bottom left panel in Fig.[TT|demonstrates that a large 
amount of scatter in excitation diagrams may be explained by 
sub-thermal excitation. The excitation temperature could there- 
fore be the kinetic temperature of dense (> 10" cm"-*) gas 
or could result from subthermal excitation of warmer gas with 
lower density. Line opacities also increase the scatter in observed 
fluxes. 

We calculate synthetic H2O spectra from RADEX'* mod- 
els of a plane-parallel slab ( van derTaket al., 2007) character- 
ized by a single temperature T, density n(H2), and H2O col- 
umn density A^(H20) with an emitting surface area A. RADEX 
is a radiative transfer code that simultaneously calculates non- 
LTE level populations and line optical depths for a plane-parallel 
slab to produce line fluxes. A large grid was calculated using 
molecular data o btained from LAMDA ( iSchoier et al.L 120051: 
lFaureetalil200l . Since this molecular data file lacks the most 
highly excited lines detected with Spitzer, individual RADEX 
models were calculated at specific gridpoints using a much 
larger and more complete database with energy levels obtained 



^ Throughout the paper all logarithms are base 10, all units for col- 
umn density are in cm"-, and all units for density are cm"^. The units 
in excitation diagrams are in number of detected molecules rather than 
column density. 

^ The terminology "warm" and "cool" components, as defined here, 
corresp onds to "hot" and "w arm" gas, respectively, in other works, in- 
cluding IVisser et"^ ( 120 llh . which considers colder gas, usually ob- 
served in the sub-mm, than the gas studied here. 



from Faure et al. (2008). 

The RADEX models are calculated to obtain a rough idea 
of the physical properties of the emitting gas. Radiative pump- 
ing is not included in the model but is likely important, espe- 
cially at low densities. Including radiative pumping would re- 
quire detailed physical and chemical modeling of the envelope 
and is beyond the scope of this work. The line profile is assumed 
to be a Gaussian profile with a FWHM of 25 kms"', based 
on the FWHM of low-excitation H2O lines observed with HIFI 
dKristensen et alj,l2010l) .^ The extinction to the warm H2O gas is 
highly uncertain and is mostly ignored (see §3.3 and Appendix D 
for a discussion of extinction estimates and their implications). 
The ortho-to-para ratio is assumed to be 3, based on flux ratios 
of optically-thin lines that range from 2.8-3.5. 

A fit (upper left panel of Fig. [T2l i to the measured fluxes 
with PACS and IRS lines with upper energy level above 400 
K (the warm component) yields acceptable solutions with high 
temperature {T > 1000 K) and H2 densities of logn < 7.5. 
Appendix C provides a detailed description of the line ratios and 
the non-detections of H^^^O emission that constrain the best-fit 
parameters. The size of the emission region further limits the 
range of acceptable parameter space (Fig.fTSli. The emitting sur- 
face area A is equivalent to the area of a circle with radius from 
25-500 AU, which is reasonably close to the projected length 
of the outflow on the sky (~ 1000 AU or 0.005 pc). If we as- 
sume that A^(H20)< 10 "* A^(H2), then the given column density 
and H2 density yields the length scales (depth along our line of 
sight) for the outflow that range from > lO"'* pc (for the accept- 
able model with the highest H2 density and lowest H2O column 
density) to > 130 pc (for the model with the lowest H2 den- 
sity and highest H2O column density). Given the projected size 
of the outflow of 0.005 pc, a length scale greater than 0. 1 pc is 
uncomfortably large and rules out solutions with log«(H2)< 5. 
The lack of obvious vibrational excitation in spectra at 6 //m 
(Maret et al., 2009; Arnold etal., 2011) limits the kinetic tem- 
perature to < 2000 K. The number of H2O molecules scales with 
density. 

Combining these analyses, we adopt the parameters T - 
1500t2000 K, logn = 6.5 + 1.5, and logA?(H20)= 17.6 + 1.2 
over an emitting area equivalent to a circle with radius 25 - 300 
AU and length scale 0.003 pc. Figure [13] shows that the PACS 
H2O spectrum is well fit with model fluxes obtained with these 
parameters. 

The temperature and column density are both inversely cor- 
related with the density, so the acceptable parameter space is 
tighter than implied by the large error bars. The choice of line 
width scales the optical depth. A broader line width would re- 
quire the same factor increase in column density A^(H20) and 
decrease in total emitting surface area. The widths of the far- 
IR lines may differ from the optically thick low-excitation H2O 
lines analyzed by Kristensen et al. (2010). The far-IR lines are 
primarily seen from the offset outflow location, while the longer- 
wavelength lines are dominated by on-source emission and in- 
clude both red- and blue-shifted outflow lobes. Very broad lines 
(> 65 kms ') are ruled out from the widths of lines in the 
Nyquist sampled spectral maps. 



http://www.strw.leidenuniv.nl/~moldata/radex.html 
' The HIFI H2O lines are much more optically-thick and likely have 
larger emitting areas than most H2O lines in the PACS spectrum, which 
may lead to differences in line profile. 



8 



Herczeg et al.: Warm water in a Class outflow 



45 
44 
43 
42 

S '^^ 

be 
o 

- 40 

39 
38 



• 5 



HgO emission lines 
T„ = 110 K 



Herschel-PACS 

Optically thin 

2<T< 10 
Optically thiclc 



Spitzer-IRS 



■ T„ = 220 K 



500 1000 1500 2000 2500 3000 
Energy (K) 



47 



46 



45 



no 
o 

^ 44 



43 



CO emission lines 



T„ = 080 K 



T„ = 280 K 



1000 2000 3000 4000 5000 6000 7000 
Energy (K) 



*: ,.. 

. : log n(H,) = 11 



03 -4 



o 



I •• • • 



T = 1500 K 



• •• I 



• • • • 



log niK.,) = 6.5 



500 1000 1500 2000 

Energy (K) 



12 



10 



a 

^ 6 



100 



^ W07 



log N(HjO) = 16.7 



log N(H20) = 17.5 * 

log N(H20) = 18.3 
log N(H5,0) = 19.0 




1000 
Temperature (K) 



Fig. 11. Excitation diagrams, in units of total number of detected molecules N divided by degeneracy g, for H2O (upper left), CO 
(upper right), and OH (lower right) emission lines detected with Herschel-PACS (circles) and Spitzer-IRS (purple squares). For the 
H2O excitation diagram, PACS data are subdivided into lines that are optically-thin (blue), moderately optically-thick (green), and 
optically-thick (black), for the H12 model. Most of the Spitzer lines are optically-thin. Lower left: H2O excitation diagrams obtained 
from RADEX models for optically thin models with T - 1500 K at high (blue) and low (red) density. At low density, subthermal 
excitation leads to cooler measured excitation temperatures and significant scatter in the level populations. 



4.2. Comparing H2O spectra for warm, subttiermal excitation 
and cool, tiiermalized gas 

The high temperature, low density solution presented here (here- 
after HI 2, with properties obtained from the;^'^ fit listed above) 
produces sub-thermal excitation of H2O, in contrast to the high- 
density (thermalized, with logn ~ 11, log A^(H20)= ^^ - Q' ™^ 
r = 170 K), low-temperature solution from lWatson et al.l (l2007l 
hereafter W07). The W07 slab has many optically-thick lines, 
which leads to a spectrum where the mid-IR H2O lines observed 
with Spitzer are brighter than those in the PACS wavelength 
range (Fig.[T4ll. In contrast, many of the far-IR lines in the H12 
parameters are optic ally- thin, so that most H2O emission escapes 
in the PACS wavelength range. As a consequence, in principle 
the far-IR H2O emission could trace a high temperature, low 
density region (HI 2) while the mid-IR H2O emission traces a 
high density, low temperature component (W07). 

The W07 model produces optically-thick emission in the o- 
H2O 818 - 7o7 and P-H2O 808 - 7i7 lines at 63.4 yum, with a flux 
in the P-H2O 808 -7i7 line similar to the observed flux. However, 
the two lines are observed in an optically-thin flux ratio both on- 
source and at the blueshifted outflow lobe. Therefore, the W07 



model cannot explain the PACS H2O emission located at the out- 
flow position. If the extinction in W07 is reduced to Ay ~ mag, 
then the two 63.4 jum lines both become two times weaker rela- 
tive to the mid-IR H2O lines. 

A comparison between the Spitzer/IRS spectrum and the 
synthetic H2O spectrum (see Fig. [15] and a further discussion 
of line ratios in Appendix C) shows that both the W07 and 
H12 models could explain the mid-IR H2O emission alone. 
Both models are able to accurately reproduce the emission in 
most detected lines, with a few notable exceptions. The P-H2O 
835 - 726 28.9 yum line flux is well reproduced in W07 but not 
H12. However, the wavelength of this line is more consistent 
with an OH line than with the H2O line. The OH rotational di- 
agram (Fig. fTTT i shows that the line flux is also consistent with 
fluxes in other OH lines with similar excitations. The inability of 
H12 to reproduce this line flux with an H2O model is therefore 
not significant. An OH line at 24.6 //m was also misidentified as 
0-H2O 863 - 836 despite neither W07 nor H12 being able to pro- 
duce flux in the H2O line. Because the synthetic fluxes of these 
lines are faint, the misidentification of this emission as H2O can 
help to drive a best fit physical parameters to an optically thick 
solution. 
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Fig. 12. Upper left: Contours of versus temperature and H2 density for several different H2O column densities (the different 
colors), with the minimum reduced ~ 4. The contours show where ;t'^^^/4 -1.2 (solid lines) and 1.5 (dashed lines), which 
roughly indicate the acceptable parameter space. The contours shown here are calculated to fits of lines with E' > 400 K. The 
asterisks show the parameters of the sub-thermal model presented here (HI 2) and the high density envelope-disk accretion shock 
model (W07). Upper right and lower panels: contours of total number of H2O molecules (upper right), radius (AU) of a circular 
emission area for H2O emission (lower left), and length scale (log pc) of H2O emission (lower right) compared with the best fit 
contours of «(H2) and A^(H20) for T - 1500 K. The asterisks indicate the parameters adopted for this paper. In the lower right 
panel, the solid line at 0. 1 pc shows the conservative upper limit to the length scale of the emission. The parameters adopted for this 
paper are located in the lower right of the acceptable parameter space so that the length scale for the emission is small. 



In models with high A^(H20) and n(H2), including W07, the 
line blend of 0-H2O 652 - 523 and 550 - 423 at 22.4 jum and the 
P-H2O 642 - 5 15 23.2 fim line are predicted to be strong but are 
not detected. Both W07 and H12 overpredict the flux in the H2O 
21.15 fj.m line. The H 1 2 model does not significantly overpredict 
any other line in the IRS spectrum, even if the best fit H12 fluxes 
are scaled to the level of the strongest IRS lines rather than to 
the fai-IR PACS lines. 



4.3. CO Excitation 

In the CO excitation diagram, a cool (275 + 30 K) component 
dominates mid-7 lines and a warm (880 + 100 K) component 
dominates high-7 lines. The uncertainty in temperature includes 
the choice of energy levels to fit for the high and cool compo- 
nents. The two excitation components could relate to regions 
with different kinetic temperatures or with different densities. 
RADEX model s of CO were run using molecular data obtained 



from LAMDA (ISchoier et all 120051: lYang et al L l2010h wifli a n 
extrapolation of collision rates up to 7 = 80 bv (lNeufeldll2012l) . 



This analysis demonstrates that a high temperature, low den- 
sity model (HI 2) of the H2O emitting region can reasonably re- 
produce both the PACS and IRS spectra. On the other hand, a 
low temperature, high density model (W07) cannot reproduce 
the PACS spectrum, is inconsistent with the spatial distribution 
of the H2O 63.4 /urn line emission, and overpredicts the emission 
in several lines in the IRS spectrum. 



The detection of high-7 CO lines with an excitation temper- 
ature of ~ 880 K requires logn(H2)> 6 for reasonable kinetic 
temperatures (< 4000 K). For T = 1500 K, logn(H2)= 6.5, and 
an emitting area with radius 100 AU, the physical parameters 
adopted to explain the water emission, produces an excitation 
temperature of 950 K for lines with J - 30 - 45. This tem- 
perature is sensitive to the density, with logn(H2)= 6.0 leading 
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Fig. 13. Segments of the PACS spectrum (solid black line). The H2O line fluxes (shaded in blue) are obtained from RADEX, CO 
line fluxes (red) from the thermal distribution in Fig. 11, and OH (green) and [O I] line fluxes (purple) from Gaussian fits. 



to an excitation temperature of 640 K, The total number of CO 
molecules, N, for log n{H2)- 6.5 is log N = 48.4. 

In principle the two temperature components could relate to 
a single region with high temperature (~ 4000 K) and log n(H2)< 
4 dNeufeld, 2012). In this case, the CO emission would be phys- 
ically unrelated to the highly excited H2O emission, and the CO 
abundance in the highly excited H2O emission region would be 
much smaller than that measured here. 



4.4. OH Excitation 

The OH excitation diagram shows cool (60 + 15 K) and warm 
(425 + 1 00 K) components . RADEX models of the low excitation 
levels of OH ( £„„ < 1000 K) wer e run, using collisional rate co- 
eflicients from (Offer et aL , 19941) an d energy levels and Einstein 



A values from (Pickett ~t all 1998h . Molecular data for higher 



4.5. O Excitation 

From the 63.18 /xm line flux and assuming T - 1500 K, the 
total number of neutral O atoms log AT - 47.9. This number 
is robust to changes in temperature and density within the pa- 
rameter space discussed here, based on RADEX models of [O 
I] lines. If the [O I] emission is spread out over a circle of 100 
AU in radius, the column density log N(0)= 16.5 is much less 
than the column density required for the line to become optically 
fliick (log N(0)~ 19 for a Gaussian profile with a FWHM of 25 
km s '). These physical properties are adopted from the H2O 
emission for simplicity but are likely incorrect because the loca- 
tion of [O I] emission is spatially different than the H2O emission 
(left panel of Fig.|6]l. 



5. DISCUSSION 



excitation levels were obtained from HITRAN (Rot hman et all 5.1. The Origin of Highiy -Excited H2O emission 

Prior to this work, H2O emission has been attributed to three 
different regions in or near the IRAS 4B environment: 

(I) A Compact Disk: Spectrally narrow (FWHM~ 1 km s ') 
HPO emission is produced in a (pseudo)-d isk with a ra dius ~ 25 



2009v) . A 1500 K gas with logn(H2)= 6.5 and emitting area of 
radius 100 AU roughly reproduces the emission in the cool com- 
ponent OH emission, with logN(OH)= 17.3. Whether these pa- 
rameters could also reproduce the highly excited OH emission is 
not clear. 

The RADEX model fluxes are also somewhat discrepant 
with the observed fluxes. The 24.6 /vm line flux is much lower 
than predicted. In addition, all detected OH doublets have sim- 
ilar line fluxes but the RADEX model predicts different fluxes 
in several transitions. A lower opacity, caused either by broader 
lines or a lower column density and larger emitting area, would 
alleviate some of these discrepancies. 

The OH molecule has energy levels with high critical den- 
sities and with strong far-IR transitions to low energy lev- 
els that are favorable to IR photoexcitation. Radiative pump- 
ing may therefore se verely alter the level populations (e.g. 
IWamnfler etani2010h . The IR pumping would increase the pop- 
ulations in excited levels, which may cause us to overestimate 
the total number of OH molecules for the given temperature and 
density. As with H2O, a rigorous assessment of the IR pumping 
requires a full physical model of the envelope and is beyond the 
scope of this work. 



AU around IRAS 4B (J0rgensen & van Dishoecki [201 Ol) . From 
the inferred li\^0 column density, the compact disk is optically- 
thick in most H^^O rotational transitions. The disk covers only a 
small area on the sky and therefore contributes very little emis- 
sion to the broad lines detected with HIFI and to the mid- and 
far-IR H2O emission. From the assumed T - 170 K and result- 
ing column density log A^(p-H20)~ 18.4 and b - I km s the 
p-Hj^O 331 - 2o2 138.5 fim line, from the same upper level as 
the observed Hj^O line, would have a flux of 3 x 10"^^ W m"^, 
50 times smaller than the measured line flux. Our PACS obser- 
vations are unable to probe this disk component. 

(2) Outflow emission, (a) in low-excitation lines with a 
large beam: Spectrally-broad emission in low-excitation H2O 
lines was observed in spatially-unresolved observations over a 
~ 20 - 40" beam and was attribut ed to outflows based on the 
line widths (iKristensen et all [201 Oh . The outflow is compact on 
the sky, so that both the red- and blue-shifted outflow lobes are 
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Fig. 14. The synthetic spectra for the H12 (blue) and W07 (red) 
models at a spectral resolution R - 1500. Because the W07 
models are optically-thick in many highly-excited levels, the 
strongest H2O lines peak in the mid-lR rather than the far-IR. 
The H12 model predicts some rovibrational emission at 6 yum 
because of the high temperature, with a strength that could fur- 
ther constrain the temperature and density of the highly excited 
H2O. 



located within this aperture. And (bj, in masers: spectrally nar- 
row (~ 1 km s"') and spatially compact H2O maser emission i s 
seen from dense gas in an outflow (e.g. iDesmurs et al.L l2009h . 
The relationship between the maser emission and the molecular 
outflow is unclear. 

(3) Disk-envelope shock: H2O emission in highly-excited 
mid-lR lines, which are spectrally and spatially-unresolved at 
low spectral and spatial r esolution, wa s attributed to the disk- 
envelope accretion shock (Wats on et al.L 120071) . The primary ar- 
gument for the accretion shock is that high H2 densities (> 10'" 
cm~^) are needed to populate the highly excited levels, which 
have high critical densities. Such high densities are expected in 
an envelope-disk accretion shock but are physically unrealistic 
for an outflow-envelope accretion shock because envelope den- 
sities are much lower than disk densities. 

Our primary goals in this work are to use the spatial dis- 
tribution and excitation of warm H2O emission to test (3), 
the envelope - disk a ccretion shock interpretation proposed by 
IWatson et al.l (l2007h . and to subsequently use the far-IR emission 
to probe the heating and cooling where the emission is produced. 
In the following subsections, we discuss the outflow origin of the 
highly-excited H2O emission and subsequently discuss the im- 
plications for the envelope-disk accretion shock and outflows. 



5.2. Outflow Origin of Higlily Excited H2O Emission 

The H2O emission detected with PACS is spatially off'set from 
the peak of the far-IR continuum emission to the south, the 
direction of the blueshifted outflow. Of the mapped lines, the 
H2O 8 18 - 7o7 and 808 - 7i7 lines at 63.4 fim are closest in ex- 
citation to the mid-IR Spitzer lines. The location of the emis- 
sion in both 63.4 fim H2O lines is consistent with the location 
of near-IR emission from IRAS 4B imaged with Spitzer/IRAC. 
The full PACS spectrum from 50-200 pm demonstrates that all 
other molecular lines are also offset from the location of the sub- 
mm continuum peak (Fig. |5]l. These lines are produced in the 
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Fig. 15. A comparison of H12 (red dashed line, scaled to the 
flux in the 63.4 pm lines) and W07 (yellow filled regions, scaled 
to the flux of the 35.5 pm line) model spectra to the Spitzer IRS 
spectrum. Green vertical lines mark the wavelengths of OH lines 
that are detected (solid lines) or are too weak or blended to be 
detected (dotted lines). Both models reproduce most SpitzerflRS 
lines reasonably well, with a few exceptions described in detail 
in the text. 



southern, blueshifted outflow lobe of IRAS 4B, as shown in the 
plots and ca r toon o f Figs. |6||2l Contemporaneous to our work, 
iTappe etal] (1201 2|) found that the spatial distribution of H2O 
emission in Spitzer/IRS spectra is consistent with an outflow ori- 
gin. The images of CO 24-23 and 49^8 also demonstrate that 
highly excited CO emission is produced in the blueshifted out- 
flow lobe. The CO 49-48 emission is located closer to the central 
source than the highly excited H2O emission. 

The combined HerschelfPACS and Spitzer/IRS excitation di- 
agram does not show any indication of multiple components in 
the highly-excited H2O lines. RADEX models indicate that the 
highly excited PACS lines are consistent with emission from 
a single slab of gas with T ~ 1500 K, logA^(H20)~ 17.6, 
log n ~ 6.5, and an emitting area equivalent to a circle with ra- 
dius ~ 100 AU. The same parameters reproduce the mid-IR H2O 
emission lines detected in the Spitzer/IRS spectrum of IRAS 4B. 
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These same physical conditions may also produce CO funda- 
mental emission, which could explain the bright IRAC 4.5 fim 
emission from the IRAS 4B outflow (Tappe et al., 2012). 

The on-source component of the far-IR H2O is not as well 
described than the outflow component because this component 
is faint in the short wavelength PACS lines, which are able to 
constrain the properties of the emission at the offset outflow po- 
sition. The on-source emission suffers from higher extinction, so 
the brightest lines are at longer wavelengths, have low excita- 
tion energies, and are optically thick. Our RADEX modeling is 
restricted to the higher excitation component of the H2O emis- 
sion. However, the similarity of the ratio of the two 63.4 jjm 
lines at the sub-mm continuum location and at the outflow lo- 
cation suggests similar excitations. The non-detection of Hj^O 
lines (see Appendix C) place a strict limit on the optical depth of 
the H2O Hnes. The H2O 108.1 jum line, of which 70% is located 
on-source, likely traces the same material as the HIFI spectra 
of low-excitation H2O lines (iKristensen et al.Ll201Ql) . The width 
of the HIFI emission (FWHM~24 km s~', with wings that ex- 
tend out to ~80 km s"') is consistent with an outflow origin and 
inconsistent with a slow (~2 kms"') envelope disk accretion 
shock that would be expected for infalling gas. 

In sub-mm line emission, both the red- and blueshifted out- 
flows are detected and spatially separated (J0rgensen et al. 2007; 
Yildiz et al. submitted), which indicates that the outflow is not 
aligned exactly along our line of sight to the central object. In 
previous near- and mid-IR imaging, only the blueshifted outflow 
is detected and the redshifted outflow is invisible. The extinction 
to the redshifted outflow is so high that even the far-IR emis- 
sion is obscured by the envelope. The outflow therefore cannot 
be aligned in the plane of the sky, and likely is aligned to within 
45° of our line of sight to the central star This large extinction is 
consistent with the outflow angle of ~ 15° relative to our line of 
sight, as inferred from assuming that the IRAS 4B outflow has 
the same age as the IRAS 4A outflow (Yildiz et al. submitted). 

The relationship between the outflow and the maser emis- 
sion is uncertain. The maser emission has a position angle of 
151° from IRAS4B (iPark & ChoiL l2007l iMarvel et al.L [200I 
iDes murs et al], l2009h . in co ntrast to the ~ 174° p osition an- 
gle of the molecular outflow (|j0rgensen et al.L l2007h . Unlike the 
molecular outflow, the maser emission is located in the plane of 
the sky (~ 77° relative to our line of sight), based on the ra- 
dial velocity and projected velocity on the sky. In this case, the 
outflow would have a dynamical age of only 100 yr. However, 
the high extinction to the redshifted outflow and low extinc- 
tion to the blueshifted outflow together suggest that the out- 
flow is aligned closer to our line of sight than in the plane of 
the sky. Jet precession, perhaps a result of binarity of the cen- 
tral object, has been suggested as a possible explanation for the 
difference between th e maser emission and molecular outflow 
(IDesmurs et al.Ll2009t) . 

H2O maser emission is produced in ga s with densities of 
10^ - lO** cm-^ (Kaufman & Neufeld','l996"). The density of the 
far-IR H2O emission could be as high as ~ 10^ cm"^. However, 
the maser emission is produced in very small spots within the 
outflow, while the H2O emission detected here is likely spread 
over a projected area equivalent to a circle with a radius ~ 100 
AU. In principle, the maser emission could simply be the small- 
est, densest regions within the same outflow that produces the 
far-IR H2O emission, although the maser and far-IR H2 emis- 
sion may be unrelated. 

In summary, we conclude that most of the H2O emission 
from IRAS 4B is produced in the blueshifted outflow. The highly 
excited H2O emission is located at the blueshifted outflow posi- 



tion. The redshifted outflow is likely not detected in these lines 
because any emission is obscured by the envelope. The optically 
thick lower-excitation H2O lines have an additional on-source 
compo nent and are likely similar to the spectrally broad outflow 
seen by IKristensen et afl(l20 10). Similarly, the cool components 
of OH and CO emission are also likely produced in the outflow 
rather than a quiescent envelope. 

5.3. Implications for Envelope-Disk Accretion 

The primary motivation for invoking the envelope-disk accretion 
shock to interpret the mid-IR H2O lines was the high critical den- 
sities for the mid-IR H2O lines. We have demonstrated that these 
lines may also be produced in gas with sub-thermal excitation. 
An additional component besides the outflow does not need to 
be invoked at present to explain the presence of the mid-IR H2O 
emission. 

The envelope-disk accretion shock may still have some un- 
detected contribution to the on-source H2O emission seen with 
PACS, in which case the mid-IR H2O emission could in princi- 
ple be the combination of the blues hifted outflow lobe and the 
envelope-disk accretion shock. The IWatson et al.l (j2007) inter- 
pretation of the H2O emission requires high line optical depths, 
including in the 63.4 /im H2O lines. However, the 63.4 yum H2O 
ortho/para lines are observed at the optically-thin flux ratio, with 
an upper limit that no more than 30% of the observed emission 
may be produced in optically thick gas. Since the envelope-disk 
accretion shock would only contribute to the on-source emis- 
sion, the maximum contribution to the total flux in the 63.4 jum 
lines is ~ 10%. In this case, we find that the luminosity of the 
disk-envelope acc retion shock would be 1.5 x 10"^ L©, 5% of 
that calculated bv lWatson et al.l (l2007l) ^. The accretion rate onto 
the disk would correspond with ~ 3 x 10"^ yr"', which is 
likely too low to represent the main phase of disk growth. 

An alternate and likely explanation is that any disk-envelope 
accretion shock is buried inside the envelope, with an extinc- 
tion that prevents the detection of any mid-IR emission pro- 
duced by such a shock. The on-source H2O and CO emission is 
likely produced in an outflow that is seen behind Ay ~ 700 mag. 
of extinction. In this case, whether the envelope-disk accretion 
shock exists is uncertain and the accretion rate is completely un- 
constrained. The prospects for studying envelope-disk accretion 
likely require a high-density (> 10" cm""*) tracer ( Blake et akl 
199£) observed at high spatial and spectral resolution in the sub- 
mm, where emission is not affected by extinction. 

5.4. Gas Line Cooling Budget and Abundances for IRAS 4B 

The bulk of the far-IR line emission is produced in the outflow. 
The emission is smoothly distributed within the outflow and is 
consistent with emission from two different positions, one lo- 
cated on-source and one located off-source. 

Table|3]lists the total cooling budget attributed to the molec- 
ular and atomic emission in the far-IR, as extracted from the 
on-source position and the blueshifted outflow lobe. Most of 
the cooling is in molecular lines rather than atomic lines, which 
is consistent with previous estimates of far-IR cooling from a 
larger sampl e of C lass objects that were observed with ISO 
( iNisini et all |2002|) . The abundance ratio^ of 0H/H20~ 0.2 
is consistent with the abundance ratio OH/H20> 0.03 for the 
outflow from the high-mass YSO W3 IRS 5 CWampfler et al.L 

* The rate is also adjusted from lWatson et alj ( i2007t) for updated dis- 
tance of 235 pc. 
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I2OI ih . the only other YSO with such a measurement. The 
H2O/CO abundance is ~ 1, on the high end of the value of 0. 1-1 
found from emis sion in lower-excitation lines of CO and H2O 
dKristensen et al.Ll2010i) . The ratio of H2O/H2 is ~ 10 based 
on the total number of H2 molecules (logA/^(H2)= 52.2) calcu- 
lated from the extinction corr ected fluxes of pure-rotational H2 
lines from (iTappe et al.Ll2012h . The H2 and H2O emission trace 
the same projected region on the sky (right panel of Fig.|6l with 
the SpitzerflRAC imaging dominated by H2 emission). 

In contrast to the luminous molecular lines, the [O I] emis- 
sion from IRAS 4B is surprisingly faint. The [O I] 63 jum line is 
typically the brightest far-lR emission line from low-mass pro- 
tostars, with an average lumi nosity of 10' -^ L/,,,/ for Class stars 
and 1 0"^ hhoi for Class I stars dGiannini et a l.. 2001; Nisini et al.L 
I2OO2 '). For IRAS 4B, the [O I] flux is 10"^^ Uoi- The weak [O 
I] emission may be a signature of the youngest protostars with 
outflows that have not yet escaped the dense envelope, as may be 
the case for IRAS 4B. Although some [O I] emission could be 
located behind the envelope for IRAS 4B and other Class ob- 
jects, the non-detection of the [O I] 145.5 fim line suggests that 
neglecting extinction does not lead to a serious underestimate 
in [O I] line lumonosity. The PACS spectrum also does not show 
any evidence for high-velocity [O I] emission from the jet, unlike 
the mo re evolved embedded object HH 46 (Ivan Kempen et all 
I2010ah . 

The abundance ratio of O/H2O is ~ 0. 1 for the detected emis- 
sion in the blueshifted outflow. The [O I] emission is predomi- 
nately located alongside the outflow axis and is offset from the 
location of the H2O emission. Within the shock that produces the 
highly excited H2O emission, the O/H2O abundance ratio must 
be much lower than 0.1. The non-detection of [C II] emission 
suggests that the ionization fraction in the shock is low. At least 
90% of the gas-phase O is in H2O, CO, and OH, whichis consis- 
tent with the high H2O/H2 ratio measured above. 

Relative to other sources, the listed H2O colum n density is ~ 
10^ - 10"* times larger th an that found from LI 157 dNisini et al.L 
I20T0I IVasta eTall I20TI . Some of this discrepancy might be at- 
tributed to methodology. The H2O column density calculated 
here is measured directly from the opacities in many different 
lines, which are spectrally unresolved. In contrast, the H2O col- 
umn density is calculated from optically thick lines, many of 
which are spectrally resolved but measured with a range of beam 
sizes. However, the H2O abundance from IRAS 4B may be much 
larger than that of most other young stellar objects, possibly be- 
cause of some physical difference in the outflow properties. 

5.5. Shock Properties for the Highly Excited Molecular 
Outflow 

The H2O line cooling is qualitatively consistent with ex- 
citation of molecular g as predicted from C-shocks by 
iKaufman & NeufeldllfToM ). The OH emission also likely traces 
shocks because OH has a high critical density and is an inter- 
mediate in the high temperature gas phase chemistry that pro- 
duces H2O. The two-temperature shape of the CO excitation 
diagram is expected fro m models of shock and UV-excitation 
of outflow cavity walls (IVisser et al.L 1201 lb and is qualitatively 



' As calculated from the highly excited H2O emission and the low 
temperature component of OH emission, both modeled with RADEX 
with the same temperature and density. It is unclear whether the warmer 
component of OH emission can be reproduced by these same parame- 
ters. The OH abundanc e may be severely affected by infrared pumping 
dWampfler et 



consistent with similar observations of the Class I sources HH 
46 and DK Cha (van Kempen et al. 2010ab). The ratio of H2O 
to CO luminosity is higher for IRAS 4B than for the Class I 
source HH 46 and for the borderline Class I/II source DK Cha. 
These differences may indicate that shock heating of the enve- 
lope plays a more im portant role than U V heating during the 
Class stage. Indeed, IVisser et"an (1201 ih suggested an evolu- 
tionary trend that for more massive and denser (younger) en- 
velopes the shock heating should dominate. The different evo- 
lutionary stages may not change the shape of the CO ladder, as 
the relative contributions of the mid-7 (15-25) and high-/ (25- 
50) CO lines are similar for DK Cha and IRAS 4B. The highly 
excited H2O emission may also be especially strong from IRAS 
4B because the outflow covers only a very compact area when 
projected on the sky. 

Highly excited H2O and OH lines were previously detected 
in mid-IR Spitzer/IRS spectra of the HH 211 bow shock, which 
is strong enough (~ 200 km s"') to dissociate molecules at the 
location of direct impact. Such a high velocity shock also pro- 
duces strong UV radiation (see models by Neufeld & Dalgarno 
1989 and observations by, e.g., Raymond et al. 1997 and Walter 
et al. 2003) that photodissociates molecules both upstream and 
downstream of the shock. In Fig. 2 of Tappe et al. (2008), many 
mid-IR OH emission lines are stronger than the H2O lines from 
HH 211. In contrast, the mid-IR H2O lines are much stronger 
than the mid-IR OH lines from IRAS 4B. In addition, the OH 
emission from HH 21 1 is detected from levels with much higher 
excitation energies than detected here, in a pattern that is con- 
sistent with prompt emission following production of OH in ex- 
cited levels through photodissociation of H2O by far-UV radia- 
tion (presumably Lyman a, Tappe et al., 2008). The atomic fine- 
structure lines are also brighter from the HH 21 1 outflow than 
from IRAS 4B dGiannini et al.L 120011: iTat^ne et al.Ll2008l) . Thus, 
UV radiation likely plays a large role in both the chemistry and 
excitation of OH at the HH 211 shock position, which is well 
separated from the YSO itself and is outside the densest part of 
the envelope. 

In the case of IRAS 4B, the combination of bright H2O emis- 
sion and faint [O I] emission suggests that the H2O dissociation 
rate is low. Moreover, the OH emission from IRAS 4B is not 
seen from the very highly excited levels detected from HH 211. 
Thus, OH likely forms through the traditional route of high tem- 
perature (> 230 K) chemistry of reactions between O and H2 to 
form OH. The OH can then collide with H2 to form H2O. These 
reactions con t rol the oxygen chemistry in dense C-type shocks 
dDraine et al.L 119831; iKaufman & Neufeldl Il996l) . The balance 
between O, OH and H2O in well-shielded regions depends pri- 
marily on the H/H2 ratio of the gas. In dense photo-dissociation 
regions, these same high temperature reactions are effective in 
forming H?0, but the strong UV field can drive H2O back into O 
and H2O dSternberg etal.Lfl99l . Indeed, for the Orion Bar, the 
very strong UV field leads to a higher OH/H2O abundance ra- 
tio (OH/H20>l) than detected here (0H/H20~ 0.1), in addition 
to emission in other diagnostics of strong UV radiation (strong 
C 111 and CH^ emission ) that are detected from the Orion Bar 
Goicoechea et al.L 12011 ') but are not detected from IRAS 4B. 
Thus, the C-shock that produces the OH and H2O emission from 
IRAS 4B is likely not irradiated by UV emission. The shocked 
gas may be shielded from any UV emission produced by the 
central source and internal shocks within the jet. 

This shielding lends support to a C-type shock explanation 
rather than molecular formation downstream of a dissociative 
J-type shock because such a shock would produce UV emis- 
sion (see discussion above). A non-dissociative shock is also 
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consistent wi t h the presence of H2 em ission from IRAS 4B 
(lArnold et all l201lHTappe et all l2012h . This scenario is dif- 
ferent than that postulated for HH 46 by (j^n Kemper Tet all 
l201Qah . where the on-source O and OH emission was thought 
to be produced by a fast dissociative J-type shock based on the 
different spatial extents of O H and H20. This scena rio also dif- 
fers from the interpretation ofWampfl er et alJ (1201 ih that the OH 
and H2O emission from the high-mass YSO W3 IRS 5 arises ei- 
ther in a J-type shock or in a a UV-irradiated C-shock. While the 
J-type shock or UV-irradiated C-type shock is an unlikely expla- 
nation for the H2O emission, the OH emission could be produced 
in a different location than the H2O emission, as is the case for 
[OI]. 



6. CONCLUSIONS 

We have analyzed //eric/ieZ/PACS spectral images of far-IR H2O 
emission from the prototypical Class YSO IRAS 4B. Table 
5 summarizes the different components of H2O emission from 
IRAS 4B. We obtained the following results: 

1. ) A rich forest of highly-excited H2O, OH, and CO emis- 
sion lines is detected in the blueshifted outflow from IRAS 4B. 
The spectrum is more line-rich than any other low-mass YSO 
that has been previously published. 

2. ) Nyquist-sampled spectral maps place the highly-excited 
63.4 jum lines at a average distance of 5'.'2 (projected distance of 
1130 AU) south of IRAS 4B. The lower-excitation H2O 108.1 
fim line has a centroid closer to the peak of the sub-mm contin- 
uum emission and has a larger spatial extent than the 63.4 fim 
emission along the outflow axis. The far-IR H2O emission can 
be interpreted as one component located at the blueshifted out- 
flow position and a second component at the peak position in the 
mass distribution (sub-mm continuum peak). The redshifted out- 
flow lobe is not detected in highly-excited H2O emission, likely 
because of a high extinction (> 1500 mag.) through to the back 
side of the envelope. 

3. ) The highest excitation lines detected with RA^CS are 
optically-thin, indicating an ortho-to-para ratio of ~ 3. RADEX 
models of the highly excited H2O lines indicate that the emis- 
sion is produced in gas described by T ~ 1500 K, log n(H2)~ 6.5 
cm""*, and logA^(H20)~ 17.6 cm"^ over an emitting area equiv- 
alent to a circle with radius ~ 100 AU. These same physical 
parameters can reproduce the mid-IR H2O emission seen with 
SpitzerflRS and the CO and OH emission seen with PACS. The 
total mass of warm H2O in the IRAS 4B outflow is about 140 
times the amount of water on Earth. 

4. ) From results (2) and (3), we conclude that the bulk of the 
far-IR H2O emission is produced in outflows. Any contribution 
of the envelope-disk accretion shock to highly-excited H2O lines 
is minimal. Moreover, at present the mid-IR H2O emission does 
not offer any support for the presence of an envelop e-disk accre- 
tion sh ock, as had been previously suggested by Watson et alj 
(l2007h . Any mid-IR H2O emission produced by the envelope- 
disk accretion shock is likely deeply embedded in the envelope. 

5. ) In the blueshifted outflow lobe over 90% of the gas-phase 
O is in H2O, CO, and OH, with H2O twice as abundant than CO 
and 10 times more abundant than OH. The cooling budget for gas 
in the envelope around IRAS 4B is dominated by H2O emission. 

6. ) The H2O emission traces high densities in non- 
dissociative C-shocks. In contrast, much of the heating of lower- 
mass envelopes occurs by energetic radiation. The highly excited 
H2O in the shock-heated gas is likely shielded from UV radia- 
tion produced by both the central star and the bow shock. The 



OH likely forms through reactions between O and H2 and pro- 
vides a pathway to form H2O. 
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Table 1. Location of PACS line and continuum emission from the IRAS 4B system 



Species 


Line 




,1 ijjm) 


ARA (") 


ADec (") 


a(RA) (") 


(r(Dec) (") 


CO 


49^8 


6457 


53.9 


3.4 ± 1.5 


-2.4 + 0.5 


1 + 1 


9.9 ± 1.0 


0-H2O 


532 - 5o5 


732 


54.5 


1.0 ± 1.0 


-5.4 ± 0.5 


9.9 ± 1.0 


10.4 ± 1.0 


[O 1] 


'Pi Pi 


228 


63.2 


-1.3 + 0.3 


-3.9 ±0.2 


11.8 ±0.7 


12.2 ±0.8 


0-H20 


§18 ~ 7o7 


1070 


63.3 


0.2 ± 0.2 


-5.1 ±0.3 


11.2±0.3 


12.2 ± 0.4 


0-H20 


808 - 7 17 


1070 


63.5 


0.2 ± 0.3 


-5.2 ± 0.4 


11.3 ±0.6 


14.4 ± 1.0 


Continuum 




63.3 


0.4 ± 0.4 


-0.7 ± 0.5 


11.5 ±0.5 


10.1 ±0.8 


0-H2O 


221 - lio 


194 


108.1 


0.3 + 0.3 


-3.1 ±0.2 


12.5 ± 0.4 


18.1 ±0.4 


CO 


24-23 


1524 


108.8 


0.0 ± 0.3 


-3.7 ± 0.3 


12.2 ± 0.5 


17.0 ±0.6 


Continuum 




108.5 


-0.05 ± 0.5 


-0.8 ±0.1 


12.5 ± 0.2 


11.3 + 0.4 


Continuum 




189.9 


(0.0) 


(0.0) 


16.5 ±0.5 


14.4 + 0.3 



Table 2. OH lines detected in the SpitzerfIRS spectrum of IRAS 46" 



Line ID 


E,„ (K) 


logA,,, (s-') 




Flux" 


Err* 


'ni/2 n3/2 7 = 9/2- 7/2^ 


875 


-1.40 


24.62 


1.6 


0.2 


2n3/2 / = 21/2- - 19/2+ 


2905 


1.29 


27.39 


6.0 


0.4 


2n3/2 7 = 21/2+-19/2- 


2899 


1.29 


27.45 


8.7 


0.4 


2ni/2 J = 19/2 - 17/2^^ 


2957 


1.28 


27.67 


8.3 


0.2 


2ni/2 -2n3/2/ = 7/2-5/2^ 


617 


-1.50 


28.94 


8.3 


0.2 


2n3/2 / = 19/2+ - 17/2- 


2381 


1.16 


30.28 


11.9 


0.3 


2n3/2 / = 19/2- - 17/2+ 


2375 


1.15 


30.35 


13.1 


0.3 


2n,/2 J = 17/2+ - 15/2- 


2439 


1.14 


30.66 


7.2 


0.6 


2ni/2 J = 17/2- - 15/2+ 


2436 


1.14 


30.71 


10.9 


0.6 


2n3/2 J = 17/2- - 15/2+ 


1905 


1.01 


33.86 


11.3 


0.8 


2n3/2 J = 17/2+ - 15/2- 


1901 


1.00 


33.95 


8.4 


0.8 


2ni/2 J = 15/2 - 13/2^ 


1969 


0.98 


34.61 


9.3 


1.0 



"Listed OH lines are the sum of unresolved triplet hyperfine structure transitions. 
* 10-^^ W cm-^, with 1 - cr error bars. 

"Includes two sets of unresolved triplets with different parities. 



Table 3. Cooling Budget with Herschel/PACS" 



Species 


Obs. 


Fraction of 


Cool Component 


Warm Component 


RADEX Model 




logL* 


Gas Cooling* 


T(K) 


logioA^" logL 


T(K) 


logio^" 


logL 




logL 




H2O 


-1.6 


0.45 


110 


47.0 -1.1 


220 


46.0 


-1.6 


(220) 


48.2 -1.5 




CO 


-1.6 


0.45 


275 


49.6 -1.8 


880 


48.5 


-1.9 


(880) 


48.5 -2.0 




OH 


-2.3 


0.09 


60 


47.1 -1.6 


425 


44.2 


-3.2 


(60) 


47.6 -2.4« 




[oiy 


-3.5 


0.005 














47.9 -3.5 





"All results from fits to excitation diagram except the luminosity of warm H2O emission. 



The warm and cool components are not necessarily the same for each molecule. 

All luminosities in units oILq. 
*PACS lines from 53-200 yum, uncorrected for extinction 
"^RADEX model with temperature of 1500 K and logn(H2)= 6.5. 
''Total number of emitting molecules that are detected. 

"The component explained by RADEX model, not the excitation temperature calculated from RADEX model. 
■'^[O I] emission is detected from a difl'erent physical location than the molecular emission. 



Table 4. Summary of components for H2O emission from IRAS 4B 



Emission source 


Probe 


Line Profile 


Opacity 


Location 


Compact Disk 


sub-mm H^^O 


Narrow, symmetric 


optically thick 


on-source 


OfF-source outflow 


mid/far-lR H2O 


Unresolved 


optically thin 


ofF-source 


Maser emission 


mm H2O 


Narrow 


high density 


off-source 


On-source outflow 


far-IR H2O 


Symmetric 


optically thin 


on-source 


Disk-envelope shock? 


(sub)-mm? 


Narrow (if present) 


optically thick 


on-source 
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Appendix A: Calibration of first and second order 
light longward of 190 fim 

PACS spectra are poorly calibrated in regions where light from 
different orders are recorded at the same physical location on 
the detector, especially between 97-103 fim and longward of 
190 jum. Most photons between 97-105 jum get dispersed into 
the second order and contaminate the flux at > 190 jum because 
of a mismatch between the grating and the filter transmission. 
Similarly, the light at 97-103 jum can be contaminated by third- 
order emission at ~ 69 /im. As a consequence, the first order 
light at 97-103 //m has low S/N and that light and the Ught at 
> 190 fj.m has not previously been flux calibrated. In this ap- 
pendix, we use PACS spectra of Serpens SMM 1 (Goicoechea et 
al., in prep.) and HD 100546 ( S turm et all 12010) . both reduced 
in the same method as IRAS 4B, to calibrate the first and sec- 
ond order emission from PACS at A > 190 //m. The continuum 
emission from HD 100546 is produced by a disk and peaks (in 
Jy) at ~ 60 /urn. The continuum emission from Serpens SMM 1 
is produced in an envelope and peaks at ~ 150 yum. 

First and second order light can be separated because the 
diffraction-limited point spread function is twice as large at 




: NGC 1333 IRAS 46; 

0.8 - 




50 100 150 200 

Wavelenth (^m) 



Fig.A.l. The ratio of flux in the central spaxel to the flux in 
the central 3x3 spaxels for HD 100546, Serpens SMM 1, and 
NGC 1333 IRAS 4B. The ratio falls linearly above > 100 yum 
until 190 jjm. Second order light contaminates the spectrum at 
> 190. Because the first and second order light have different 
point-spread functions, the plotted ratio determines the fraction 
of first and second order photons versus wavelength. For NGC 
1333 IRAS 4B, the continuum flux in the central spaxel falls at 
< 70 fim, likely because of extinction. 



200 fj.m as at 100 jum. Figure [ATI shows the fraction of flux in 
the central spaxel divided by the flux in the central 3x3 spax- 
els. This fraction decreases smoothly at > 100 jim. The fraction 
starts to rise at ~ 190 fim because of a contribution from second 
order emission. At each wavelength, this fraction directly leads 
to a ratio of first order photons to second order photons recorded 
by PACS. The second order light is then calibrated by measuring 
the first order light between 90-1 lOyum for HD 100546. Most of 
the light recorded at > 190 fim from HD 100546 is second order 
light. 

This calibration is then applied in reverse to the spectrum of 
Serpens SMM 1 to subtract the second order emission, leaving 
only first order emission at > 190 yum. A SPIRE spectrum of 
Serpens SMM 1 (Goicoechea et al., in prep) is then used to flux 
calibrate the PACS spectrum from 190-210 fim. Despite the red 
spectrum of Serpens SMM 1, ~ 50% of the photons recorded at 
202 yum are second-order 101 yum photons. 

To calibrate the PACS spectrum of IRAS 4B, the first and 
second order light are initially separated based on the point 
spread function for both wavelengths (Fig. lA.lb . The separate 
first and second order spectra are then flux calibrated using the 
relationships calculated above. Lines are identified as first or sec- 
ond order emission after searching for the correct line identifica- 
tion at A and A/2. 

Figure lA!2l shows the resulting spectrum of IRAS 4B at 95- 
105 fim. Analysis of a large sample of PACS spectra from the 
DIGIT program (P.I. N. Evans) suggest that our flux calibration 
has an uncertainty of ~ 20% between 98-103 fim, in addition 
to other sources of uncertainty in the standard PACS flux cal- 
ibration. The line fluxes measured from second order light at 
> 190 yum are consistent to 20% of the fluxes measured directly 
at 97-103 yum but with smaller error bars. The accuracy of the 
flux calibration at > 190 fim has not been evaluated but is likely 
uncertain to ~ 40%. 

The best place to observe lines between 98-103 yum with 
PACS is in the second order, despite overlap with first-order light 
because of sensitivity and higher spectral resolution. The PACS 
sensitivity to first order photons at > 200 fim is low. This new 
calibration allows us to improve the accuracy of flux measure- 
ments for lines between 98-103 yum and at > 190 fim. 

Appendix B: Unidentified Lines 

All strong lines in the IRAS 4B spectrum are identified. Table 
IB.ll lists several possible lines that are detected at the ~ 3cr sig- 
nificance level. The position of several of these tentative detec- 
tions do not correspond to expected emission lines. Although 
no o-H^^O are clearly detected, flie o-H^^^O In - lio 109.346 
fim line is expected to be among the strongest H2'^0 lines and 
is tentatively detected. The centroid of the detected emission is 
-170 km s"' from the expected centroid, based on the measured 
wavelengths of other nearby lines. Emission is detected consis- 
tent with the location of the P-H2O 551 - 624 71.787 fim line, but 
the line flux is expected to be very weak. 

These lines are all very weak and may be statistical fluctua- 
tions in the spectrum rather than significant detections of uniden- 
tified lines. 

Appendix C: Exploring the parameter space in 
model fits to H2O line fluxes 

In §4, we characterized the physical properties of the H2O emis- 
sion region by comparing the observed line fluxes to synthetic 
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100 1000 100 1000 4 6 B 10 12 

Temperature (K) Temperature (K) log n(H2) (em"'') 



Fig. C.l. H2O fluxes and acceptable parameter space of temperature and H2 density. The upper left panel shows that in the H12 
model, all far-IR H2O lines expected to be strong are detected (red dots), while those expected to be weak are undetected (black 
dots). Several lines at ~ 100 //m were likely undetected because the S/N is degraded in that wavelength region. The upper middle 
panel shows the acceptable contours for flux ratios of 0-H2O Sig - l^nH^M - 615 (blue), 0-H2O 7o7 - 616/615 - 5o5 (red), and both 
(purple), for two different column densities. The acceptable parameter space shows the Icr error bars with a 20% relative flux 
calibration uncertainty. The lower left panel shows that the non-detection of the 0-H2O Sso - 54i and 652 - 643 lines relative to 0-H2O 
616 - 5o5 line places a limit on the optical depth of the slab. The lower middle panel shows where the flux ratio of 0-H2O 845 - 734 
35.669 jum to 0-H2O 7o7 - 616 71.947 yum is equal to the observed value of 0.07. The lower right panel shows the parameter space 
(high density and high column density) ruled out by the non-detection of the o-Hj^^O 32i - 2i2 line, for four different temperatures 



fluxes calculated in RADEX models of a plane-parallel slab. In 
this appendix, we describe how the measured line fluxes con- 
strain the parameters of the highly excited gas. 



Table B.l. Unidentified lines" 




100 102 

Wavelength (/xm) 

Fig. A.2. The first and second order spectrum of IRAS 4B at ~ 
100 jum. The second order light at 100 yum (blue), calibrated from 
the ~ 200 yum spectrum, matches the standard first and second 
order spectra (red and green). The lines in the blue spectrum are 
artificially weak, by a factor equivalent to the ratio of first to 
second order light, because the lines and continuum are treated 
separately. The dotted vertical line shows the location of the CO 
13-12 line at 200.27 yum. 





Flux'^ 


Error* 


Possible ID 




A,„, (jum) 


57.930 


0.38 


0.10 










60.863 


0.27 


0.14 










71.809 


0.34 


0.09 


P-H2O 


551 - 


624 


71.787 


11.129 


0.23 


0.08 


0-H2O 


752 - 


743 


77.761 


95.799 


0.54 


0.13 










109.301 


0.39 


0.13 


o-HfO 


2,2- 


lio 


109.346 


112.846 


0.24 


0.12 


0-H2O 


441 - 


5|4 


112.802 


127.164 


0.21 


0.08 










129.885 


0.25 


0.08 










147.380 


0.28 


0.08 










170.017 


0.26 


0.07 


P-H2O 


633- 


624 


170.138 


181.595 


0.32 


0.14 










"These lines may not be real. 

*10-2° Wcm-2 s-'. 



The backbone 0-H2O lines 9()9 - 818, Sig - 1qi, 1qi - 616, 
616 - 5o5, 5()5 - 4i4 and their P-H2O counterparts are amongst 
the best diagnostics for the excitation and optical depth of the 
H2O emitting region. The relative flux calibration between these 
ortho/para pairs should be accurate to better than 5% because 
the lines are located very near each other. The flux ratio ranges 
from 2.8-3.5 in the three of the backbone line ratios with high- 
est excitation, indicating that the ortho-to-para ratio of ~ 3 is 
thermalized and that these lines are optically-thin. The 616 - 5o5, 
5o5 - 4i4, and 4i4 - 3o3 and their P-H2O counterparts have flux 
ratios of ~ 2.3 - 2.5, suggesting that either the modeled emis- 
sion is moderately optically thick in these lines, or that a second, 
optically-thick component contributes some emission to the flux 
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in lower-excitation, longer wavelength lines. The top panels of 
Figure IC.ll show the acceptable space of n(H2) and T for two 
different values of A^(H20) based on two of these line ratios. 

The statistic automates this type of analysis over all of 
the applicable lines to find the best-fit parameters for a sin- 
gle isothermal slab. Limiting the calculation to lines with 
E' > 600 K yields two acceptable parameter spaces, one at the 
H12 location (r ~ 1500 K, logn(H2)~ 6.5, and log A^(H20)~ 
17.6) and one at low-temperature, high density, and low column 
density (hereafter Xll, with T ~ 200 K, log«(H2)~ 11, and 
log A^(H20)~ 14.0). The Xll parameters differ from W07 be- 
cause W07 includes a large H2O column density to produce op- 
tically thick lines. The XI 1 parameters are disqualified because 
several lines, such as 0-H2O 330 - 221 66.44 fim, have a syn- 
thetic flux that is much stronger than the observed emission.** 
As a consequence, the x^ statistic for all lines with E' > 400 
K yields acceptable line fluxes only around the H12 parameters, 
adopted in this paper, and disqualifies the Xll parameter space. 

The 0-H2O transitions 550 - 54i at 75.91 fim and 652 - 643 at 
75.83/vmhavea2cr upper limit on the combined flux of 3x10"^' 
W cm"^. Although these upper limits are not included in the 
X^ calculations, the non-detections confirm that the parameter 
space with high H2 density and high H2O column density cannot 
explain the bulk of the far-IR H2O emission (lower left panel of 
Fig. IC.ll ). Similarly, the line ratio 0-H2O 845 - 734 35.67 fim to 
0-H2O 7o7-6i6 71.95 yields physical parameters consistent with 
the best-fit solution (lower right panel of Fig. IC.ll ). 

The non-detection of emission in H'^O lines (with the pos- 
sible exception of o-Hj^O 2i2 - lio 109.35 fim) also places a 
limit on the opacity of Hj^O lines. The bottom right panel of 
Fig. ICU shows where the flux ratio of o-H^^^O 32i - 2i2 75.87 
fj.m to 0-H2O 616 - 5o5 82.03 fj.m becomes > 50, large enough 
that the Hj^O line would be detected. This upper limit rules out 
the parameter space of high density and high column density. 

I n the H12 model, the strength of the 6 jum continuum emis- 
sion ('M aret et all '2009) is similar to the strength of the 6 jim 
rovibrational H2O lines. The non-detection with R = 60 spectra 
is marginally consistent with the predicted emission and rules 
out temperatures higher than ~ 2000 K. An Ay = 20 mag. would 
reduce the predicted 6 yum emission by 35%. Decreasing the tem- 
perature from 1500 K to 1000 K would reduce the predicted 
6 fim emission by a factor of 4.8. Measurements of emission 
in H2O vibrational lines would place significant additional con- 
straints on the properties of the H2O emission. 



Applying an extinction correction under these assumptions 
increases the excitation temperature of the warm (cool) CO com- 
ponent from 850 to 950 K (250 to 280 K). If all the CO emission 
were located behind the Av=700 mag., then the temperature dif- 
ference would become much larger, but the short wavelength CO 
lines would be very faint. The H2O temperature does not change 
significantly because the fits are based mostly on lines at short 
wavelengths, so the emission from the heavily embedded region 
is mostly extinguished. An Ay - 100 mag. to the blue outflow 
lobe would increase the relative luminosity of lines at 25 //m 
compared with those at 100 yum by a factor of four, thereby in- 
creasing the warm H2O excitation temperature from 220 K to 
240 K. An Ay > 200 mag. to the outflow is ruled out by the 
significant increase in scatter in the excitation diagram. 



Appendix D: Effect of Extinction Estimates on 
Rotational Diagrams 

Extinction can affect the derivation of excitation temperatures 
and line luminosities, especially at shorter wavelengths where 
lines typically have upper levels with higher energies. In this 
analysis, we assume that for all H2O lines, 70% of the luminos- 
ity is produced on source behind Ay - 700 mag. and 30% is 
produced off-source at the blueshifted outflow lobe, with Ay = 
mag. (see also §4.3). This description is consistent with the 
different observed spatial distributions of the H2O 108.1 and 
63.4 fim lines. The excitation conditions and emission line ra- 
tios may instead vary smoothly with distance in the outflow. 



The opposite case, where the synthetic emission is weaker than 
observed in a low-excitation line, would not disqualify a model be- 
cause low-excitation line emission may include an additional contri- 
bution from cooler gas. 
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Fig.D.l. The continuum-subtracted PACS spectrum of NGC 1333 IRAS 4B from 55-100 jum. The marks identify lines of H2O 
(blue), CO (red), OH (green), and [O I] (purple). 
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Fig.D.2. The continuum-subti-acted PACS spectrum of NGC 1333 IRAS 4B from 103-195 yum. The marks identify Hnes of H2O 
(blue), CO (red), and OH (green). 
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Table D.l. Detected lines in Herschel/PACS spectrum of IRAS 4B 



Species 


Line 


£' (K) 


logA,,, (s-'j 


(/'111) 


-lobs (/"111) 


Mux" 


Lrror" 


CO 


49-48 


6724 


-2.17 


53.897 


53.893 


0.66^ 


0.11 


0-H20 


532 - 5o5 


732 


-1.42 


54.506 


54.495 


0.61 


0.12 


CO 


48-47 


6457 


-2.18 


54.986 


54.965 


0.39 


0.07 


0-H20 


827 ~Tl6 


1274 


0.28 


55.131 


55.146 


0.53* 


0.13 


P-H20 


431 - 322 


552 


0.16 


56.325 


56.344 


0.59 


0.08 


P-H20 


9|9 - 808 


1324 


0.40 


56.771 


56.786 


0.21* 


0.07 


0-H20 


9o9 - 8 18 


1323 


0.39 


56.816 


56.822 


0.74 


0.09 


CO 


46-45 


5939 


-2.22 


57.308 


57.317 


0.27* 


0.09 


P-H20 


422 - 3l3 


454 


-0.42 


57.636 


57.655 


0.97 


0.11 


P-H20 


817 ~ ^26 


1270 


0.21 


57.709 


57.730 


0.43 


0.08 


CO 


45-44 


5688 


-2.23 


58.547 


58.570 


0.39 


0.09 


0-H20 


432 - 321 


550 


0.14 


58.699 


58.721 


1.60 


0.09 


CO 


44-43 


5442 


-2.25 


59.843 


59.851 


0.70 


0.08 


P-H20 


726 - 6|5 


1021 


0.13 


59.987 


60.003 


0.66 


0.07 


CO 


43-42 


4202 


-2.27 


61.201 


61.215 


0.79 


0.10 


P-H20 


431 - 4o4 


552 


-1.61 


61.808 


61.838 


0.40 


0.18 


CO 


42-41 


4967 


-2.29 


62.624 


62.644 


0.80* 


0.23 


0-H20 


9l8-9o9 


1552 


-0.38 


62.928 


62.949 


0.94 


0.18 


[0 I] 


3P,-3P2 


228 


-4.05 


63.185 


63.193 


1.83 


0.08 


0-H20 


8 18 - 7o7 


1070 


0.24 


63.323 


63.333 


1.92 


0.14 


P-H20 


808 - 7|7 


1070 


0.24 


63.457 


63.480 


0.68 


0.13 


CO 


41-40 


4737 


-2.31 


64.117 


64.127 


1.08 


0.13 


OH 


2n3/2 7 = 9/2--7/2+ 


512 


0.11 


65.131 


65.145 


0.99 


0.12 


0-H20 


625 - 5 14 


795 


-0.03 


65.166 


65.180 


2.94 


0.13 


OH 


2n3/2 7 = 9/2+-7/2- 


510 


0.10 


65.278 


65.297 


1.05 


0.12 


CO 


40-39 


4513 


-2.34 


65.686 


65.689 


0.76 


0.10 


0-H20 


7l6 - 625 


1013 


-0.02 


66.092 


66.105 


1.38 


0.10 


0-H20 


330 - 221 


410 


0.09 


66.437 


66.446 


3.07 


0.13 


P-H20 


331 - 220 


410 


0.09 


67.089 


67.102 


1.44 


0.09 


0-H20 


330 - 3o3 


410 


-2.07 


67.268 


67.287 


1.24 


0.10 


CO 


39-38 


4294 


-2.36 


67.336 


67.355 


1.18 


0.10 


CO 


38 - 37 


4080 


-2.39 


69.074 


69.075 


1.57 


0.12 


0-H20 


827 ~ 818 


1274 


-0.48 


70.702 


70.707 


0.45* 


0.12 


CO 


37-36 


3872 


-2.41 


70.907 


70.916 


1.29 


0.10 


P-H20 


524 - 4l3 


598 


-0.18 


71.067 


71.072 


1.98 


0.09 


OH 


J = 7/2- - 5/2+ 


617 


0.01 


71.170 


71.172 


0.60 


0.07 


OH 


J = 7/2+ - 5/2- 


617 


0.01 


71.215 


71.217 


0.75 


0.07 


P-H2O 


7 17 - 606 


843 


0.07 


71.539 


71.542 


1.44 


0.10 


0-H2O 


7o7 - 616 


843 


0.06 


71.946 


71.961 


3.79 


0.12 


CO 


36-35 


3669 


-2.44 


72.843 


72.855 


1.53 


0.11 


CO 


35-34 


3471 


-2.47 


74.890 


74.896 


2.01 


0.11 


0-H2O 


725 - 634 


1125 


-0.59 


74.944 


74.951 


0.89 


0.11 


0-H2O 


321 - 2i2 


305 


-0.48 


75.380 


75.395 


9.35 


0.25 


0-H2O 


854 - 845 


1805 


-0.24 


75.495 


75.491 


0.33 


0.10 


CO 


34-33 


3279 


-2.50 


77.059 


77.068 


2.76 


0.12 


0-H20 


423 - 3l2 


432 


-0.32 


78.742 


78.755 


7.57 


0.23 


P-H20 


615 - 524 


781 


-0.34 


78.928 


78.936 


1.09 


0.11 


OH 


'ni/2 -- n3/2 J = 1/2- - 3/2+ 


181 


-1.44 


79.115 


79.118 


2.53 


0.12 


OH 


-2 n3/2 J = 1/2+ - 3/2- 


181 


-1.44 


79.178 


79.182 


2.16 


0.12 


CO 


33 - 32 


3092 


-2.53 


79.360 


79.370 


2.66 


0.12 


P-H2O 


726 - 7l7 


1021 


-0.61 


81.215 


81.215 


0.38 


0.07 


CO 


32-31 


2911 


-2.56 


81.806 


81.817 


2.64 


0.09 


0-H2O 


616 - 5o5 


643 


-0.13 


82.031 


82.038 


6.67 


0.16 


P-H2O 


606 - 5 15 


642 


-0.15 


83.283 


83.287 


2.75 


0.11 


OH 


^U„2 J = 7/2+ - 5/2- 


291 


-0.28 


84.420 


84.418 


5.90 


0.22 


CO 


31-30 


2735 


-2.60 


84.410 


blend with OH 2H3/2 J 


= 7/2+ - 5/2" 


OH 


2n3/2 7 = 7/2- -5/2+ 


290 


-0.28 


84.596 


84.607 


3.06 


0.11 


0-H2O 


7 16 - 7o7 


1013 


-0.67 


84.766 


84.771 


0.76 


0.09 



"10 W cm ^, extracted from two spaxels as discussed in §2. 



Listed errors are Icr and do not include calibration uncertainty. 
*nux measured from outflow spaxel only. 
''Flux measured from Nyquist-sampled map. 
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Species 


Line 


E (K) 


10gA„;(S ) 


/Ivac (Mm} 


/lobs O^m) 


Flux" 


Error^ 


CO 


30 - 29 


2565 


-2.63 


87.190 


87.190 


3.38 


0.13 


P-H2O 


322 - 2ll 


296 


-0.45 


89.988 


89.988 


4.56 


0.10 


CO 


29 - 28 


2400 


-2.67 


90.163 


90.165 


4.12 


0.10 


CO 


16 — 11 


ZZ40 


-Z. / 1 


AQ Q /I A 

93.34y 


AQ 


/i Q A 

4.36 


A 1 

0.12 


P-H2U 




/o 


-U.uZ 


A/1 OAA 


y4.zy / 


A 


A 1 A 

U.IU 


0-rl2U 


O25 - OI6 


/yj 


-U. /O 


QA A/I '2 

y4.o4j 


y4.o4j 


I.JO 


All 

U.ll 


, , II f ^ 

0-rl2^ 


A A 
^41 - ^32 


/uz 


-U.oZ 




y^f. / ID 


1 nn 


n in 
u.iu 


1 1 , f\ 


^ , . A,. . 
3 15 %4 




-U. 


AC (COf; 


95.637 


4.44 


U. IZ 


Utl 


2tt 2 tt j a /0+ c /o- 


Z/0 


-z.Uj 


A^ 1 

yo.z / 1 


96.316 


0.57 


All 
0. 1 1 


Utl 


2tt 2 tt / 1 /o- ^ /0+ 
Hi/2 - ii3/2 J - il^ - 5/2 


zoy 


-z.Uj 


AA QAO 


96.394 


0.44 


AIT 
0. 12 


LU 


2 / — ZD 


ZOoO 


-I. /J 


AA T7Q 

yo. / /3 


96.775 


4.28 


A 1 Q 

0.13 


OH 


ni/2 y = 5/2 - 3/2+ 


Ate 

415 


-0.45 


98.736 


98.728' 


2.46 


A /1 1 

0.43 


1 1 

OH 


^ni/2 / = 5/2+ - 3/2 


415 


n AC 

-0.45 


AO n/ZA 

98.764 


blend with OH ^Hwz / = 


c //>— /J /0+ 

5/2 - 3/2 


O-H2O 


5o5 - 4i4 


468 


-0.41 


99.492 


99.496" 


9.85 


0.26 


CO 


26 - 25 


1937 


-2.80 


1 AA /I ^ A 

100.460 


100.464'' 


3.86 


0.20 


0-H20 


5 14 -423 


574 


-0.81 


100.912 


100.910' 


4.56 


0.22 


P-H20 


220 - 111 


1 nc 

195 


-0.58 


1 AA A01 

100.962 


100.981' 


5.42 


0.17 


P-H20 


Ol5 - O06 


TO 1 

781 


-0.87 


103.916 


103.944 


3.50 


0.80 


CO 


zj — z4 


i 794 


o /I 
-Z.S4 


1 A/i /i /I C 

104.445 


104.464 


4.22 


A 1 

0.18 


0-H2O 




194 


-0.59 


1 AO A"71 

lOs.0/2 


108.097 


9.47 


0.19 


CO 


Z4 — z3 


iODO 


o on 
-Z.isy 


1 AO TAQ 

lOo. /03 


108.774 


4.86 


A 1 C 

U.15 


I I l~\ 

P-H20 




j9o 


-U.9z 


111 AOT 

lll.oz/ 


111.646 


0.71 


A 1 A 
U.IU 


, , I I t~\ 

0-H2O 


'43 - '34 


1339 


A 
-U.O/ 


110 CIA 

11Z.51U 


112.492 


0.21* 


A AA 
U.U9 


CU 


Zo — zz 


ljZ4 


-z.y4 


I 1 '2 /ICQ 

I I J.458 


113.470 


4.88 


A 1 C 

U.15 


0-ri20 


4l4 - 3o3 


323 


A 1 
-U.Di 


1 1 Q CQA 

1 13.5 j6 


113.549 


11.70 


A 1 C 

U.15 


P-H2U 


c c 

-^33 ~ -^24 


/Id 


A no 


1 1 Q A/1-7 

1 13.94/ 


113.937 


0.31 


A 1 C 

U.15 


CD 


ZZ — Zl 


139/ 


Q AA 

-3.UU 


110 col 
1 10.581 


118.595 


6.46 


All 
U.ll 


yjti 


H3/2 J - 5/2 - 3/2 


1 or\ 
izU 


A O/C 


1 1 A 111 

1 19.233 


119.241 


2.23 


A 1 O 

U. 12 


OH 


^n3/2 7 = 5/2+ - 3/2 


120 


-0.86 


119.441 


119.447 


2.83 


0.11 


0-H2O 


432 - 423 


550 


-0.91 


121.721 


121.740 


1.17 


0.08 


CO 


21 - 20 


1276 


-3.05 


124.193 


124.204 


7.02 


0.10 


P-H2O 


4o4 - 3 13 


319 


-0.76 


125.353 


125.383 


4.73 


0.10 


P-H2O 


331 - 322 


410 


-1.11 


126.713 


126.724 


0.46 


0.08 


0-H2O 


725 - 7l6 


1125 


-0.87 


127.883 


127.922 


0.30 


0.05 


CO 


20 - 19 


1160 


-3.11 


ion o /' r\ 

130.369 


130.393 


7.38 


0.05 


0-H20 


423 - 4l4 


432 


-1.09 


132.407 


132.437 


2.73 


0.05 


0-H20 


5 14 - 5o5 


574 


-1.12 


134.934 


134.964 


1.27 


0.07 


0-H20 


330 - 321 


410 


-1.18 


136.495 


136.520 


1.63 


0.06 


CO 


19 - 18 


1050 


-3.18 


137.196 


137.230 


7.66 


0.07 


P-H20 


3l3 - 2o2 


204 


-0.90 


138.527 


138.536 


6.96 


0.07 


P-H20 


4l3 - 322 


396 


-1.48 


144.517 


144.533 


2.04 


0.07 


CO 


18 - 17 


945 


-3.24 


144.784 


144.794 


8.13 


0.09 


P-H20 


431 - 422 


552 


-1.09 


146.923 


146.945 


0.49 


0.09 


CO 


17 — 16 


A£ 

846 


-3.32 


153.267 


153.282 


8.87 


0.06 


p-rl2D 


ill - J 13 


Zyb 


-l.zo 


1 CA 1 AQ 

156. lyj 


156.190 


2.16 


A AA 
U.U6 


0-n20 


J23 - '^32 




1 AC 

-1.95 


1 C A OA /I 

156.264 


156.261 


1.55 


A AA 
U.U6 


, , II 1 \ 

0-ri2U 


c c 

J32 - J23 


lil 


1 AA 


1 AA CAO 

16U.5Uy 


160.513 


0.66 


A AT 
U.U/ 


CO 


io — lo 


TCI 


Q QA 

-3.39 


1 AO 010 

102.612 


162.825 


9.13 


A OA 

U.2U 


f\ I I 

OH 


2tt T "3 1 /T- 

II3/2 J = 3/2 - 1/2 


z/U 


1 1 A 

-1.19 


163.123 


163.149 


1.01 


A AO 

U.U8 


OH 


^n3/2 7 = 3/2 - 1/2+ 


269 


-1.19 


163.395 


163.423 


0.88 


0.04 


P-H2O 


624 - 615 


867 


-1.12 


167.034 


167.027 


0.44 


0.06 


CO 


15 - 14 


663 


-3.47 


173.631 


173.643 


9.45 


0.11 


0-H20 


3o3 - 2l2 


196 


-1.30 


174.624 


174.625 


9.53 


0.11 


0-H20 


2l2 - loi 


114 


-1.25 


179.525 


179.532 


8.86 


0.10 


0-H20 


221 - 2l2 


194 


-1.51 


180.487 


180.491 


3.56 


0.09 


CO 


14-13 


580 


-3.56 


185.998 


186.014 


8.64 


0.11 


p-H,0 


4|3 -4,14 


396 


-1.43 


187.109 


187.102 


0.93 


0.10 


CO 


13 - 12 


503 


-3.66 


200.271 


200.262 


6.33 


0.10 



cm ^, extracted from two spaxels as discussed in §2. 
Listed errors are Icr and do not include calibration uncertainty.. 

''Flux measured from outflow spaxel only. 

'Measured from first or second order emission at > 190 \xm. 
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